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RIGHT AND LEFT JOINT SYSTEM REPRESENTATION OF 
A RATIONAL MATRIX FUNCTION IN GENERAL POSITION 

(SYSTEM REPRESENTATION THEORY FOR DUMMIES) 
V.E. Katsnelson 

For a rational k x k matrix function R of one variable in general position, the ma- 
trix functions R{x) ■ R"^{y) and R~^{x) ■ R{y) of two variables are considered. For these 
matrix functions of two variables, the representations which are analogous to the system 
representation (or realization) of a rational matrix function of one variable are constructed. 
This representation of the function R{x) ■ R~^{y) (of the function R~^{x) ■ R{y)) is said to 
be the joint right (respectively the joint left) system representation of the matrix functions 
R, R^^. In these representations there appear diagonal nxrt matrices, A-p — diag (Ai, . . . , A„) 
(called the pole matrix for R) and A^/ = diag (/ii, . . . , /i„) (called the zero matrix for 
R), where A/,...,A„ are poles of i?, ^i,...,^n are poles of i?^^; and k x n matrices 
F-p and Fj\f (called the left pole and zero semi-residual matrices) and n x k matrices 
G-p and Gm (called the right pole and zero semi-residual matrices) which can be intro- 
duced from the additive decompositions R(z) — R{oo) + F-p{zI — Ap)^^G-p, R~^{z) = 
R{oo)^^ + Fj^{zl — Aj^)^^Gj\/. The right joint system representation has the form 
rIx) ■ R-\y) = I + {x -y)Fv{xI ~ Ap)-^S'')-HyI - Aj^y^Gj^, the left one has the 
form R-\x) ■ i?(y) = I + {x - y)Fj^{xI - Aj^)-^{S')-\yI - A-p)-^Gp. The n x n matri- 
ces and 5' (the so-called right and left zero-pole coupling matrices for R) are solu- 
tions of the appropriate Sylvester-Lyapunov equations. These matrices are mutually inverse: 
S"^ • 5' = • S"- = /. 

These results are essentially not new: they could be easily derived from known results 
on realization of a rational matrix functions (for example, from results by L. Sakhnovich or 
J. Ball, I. Gohberg, L. Rodman), however the method is new, as well as the emphasis on 
"the left, the right and their relationships". The presentation is oriented to a "traditional" 
analyst. No previous knowledge in realization theory of matrix functions or its ideology is 
assumed. One of the purposes of this paper is to provide a realization theory background 
for investigations of the deformation theory of Fuchsian differential system and of rational 
solutions of the Schlesinger system. As an application we also consider the spectral (Wiener- 
Hopf) factorization. 

The concluding Section 5 contains some historical remarks highlighting the role of 
M.S. Livsic as the forefather of the system realization theory. 

NOTATIONS. 

— . — def 

• C stands for the complex plane; C is the extended complex plane: C = C U oo; 

• DJlk stands for the set of all A; x A: matrices with complex entries; 

• / stands for the unity matrix of the appropriate dimension; 

• TZ{DJlk) stands for the set of all rational 9Jlfc-valued functions R with deti?(z) ^ 0. 

• V{R) stands for the set of all poles of the function i2, Af{R) stands for the set of all poles of 
the function R'^; 'P{R) is said to be the pole set of the function R, M{R) is said to be the zero 
set of the function R. 
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0. PREFACE 



The problem which wc set as a goal in this paper for scalar (i.e. complex valued) func- 
tions means to restore a rational function from its poles and zeros. The traditional solution of 
this problem uses products constructed from the poles and zeros of the function. For rational 
functions in general position, this solution can be explained particularly clearly. 

Namely, let r be a rational function in general position -"^j with the pole set V{r) and the zero 
set J\f{r). These sets V{r) and 7V(r) do not intersect (i.e. V{r) r\J\f{r) = ) and are of the same 
cardinality: #7^(r) = ij^M{r). The function r admits the representation 

r(z)=c( n (^-w)) • I n i'-^^A , (0-1) 

where c = r(oo). This mAiltipKcative representation recovers the function r from its pole and 
zero sets and from the value r(oo). Inversely, given two finite non-intersecting sets V and J\f 
{P, A/" G C, n A/" = 0) of the same cardinality and a complex number c 7^ 0, oo, we define the 
function r by the formula (0.1). This function r is a rational function in general position, the 
given sets V and J\f are its pole and zero sets P{r) and J\f{r) and c = r(oo). 

However, in view of non-commutativity of the matricial multiplication, the multiplicative 
representation (0.1) seems to be unsuitable for generalization to matrix functions. We present 
now such a representation of a rational matrix function (in general position) from its poles and 
zeros which can be generalized to the matricial case. This is the so-called system representation 
of a rational function. 

So, let again r be a rational function in general position, with the pole and zero sets 'P(r) 
and Af{r). We derive its system representation. We assume for simplicity that the function r is 
normalized by the condition r(oo) = 1. We start from the additive decomposition of the matrix 
function r: 

The condition 

rifip) = (V//p G J\f{r)) (0.3) 
leads to the system of linear equations 

E 17^ = -^ (V/^peA/-(r)). (0.4) 

x,ev{r) ^1 

Thus, to restore a rational function in general position from its poles and zeros, we have to 
solve the linear system (0.4) with respect to and then to substitute these into (0.2). Since 

#P(r) = 4f^M{r) =^ n, the matrix S of the system (0.4) is square: 

S = ||*p,q||l<p,q<n) ^P,<1 ~ (^•^) 

/ip — Aq 

The system (0.4) is uniquely solvable: its determinant (which is known as the Cauchy determi- 
nant) can be calculated explicitly (see, for example, [PS], Pt.VII: §1, no. 3). From this explicit 
expression for the determinant it is evident that det S 0. 



^This means that all poles and zeros of the function r are simple and r(oo) 7^ 0, 00. 
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We can formulate this method of restoring of the function r from V{r) and J\f{r) in the 
matricial form. Let A-p and Aj\f be the diagonal matrices constructed from V{r) and M{r): 

= diag(Ai, A2, . . . A„), = diag(/xi, //2, ■ ■ ■ (0.6) 

Let e be the n-row (i.e. 1 x n matrix): 

e=[l, 1, (0.7) 

As usual, by e* we denote the Hermitian conjugate to e: e* is a n-column (i.e. n x 1 matrix). 
The representation (0.2) can be put down in the form^ 

r{z) = l + e{zl- Ar)-'[^i, 6, • • • , ^n]^- (0.8) 

The system (0.4) can be presented in matricial form: ^[^i, ^2, • • • , ^n]^ = —e*. Thus, 

[^1,6, ...,^nf = -S-'e*. 

Substituting this expression for [^1, ^2, • • • , ^n]'^ into (0.8), we come to the representation 

r{z) = l-e{zI-A'p)-^S-^e*. (0.9) 

In the same way we can obtain the representation for the function r~^. Starting from the additive 
representation 

r-^(^) = l+ ^ (0.10) 



we come to the linear system (with respect to rjp) 

tipeMir) 

which can be put down in matricial form: 

[m, V2, Vn]S = e, 

or 

iVi, V2, Vn] =eS'\ 

Thus, 



r-\z) = l + eS-\zI - Aj^y^e*. (0.11) 

Here the matrix S is the same as in (0.5) ! 

The matrices S, A-p, Aj^j- are involved in the equality 

Aj^S - SA-p = e*e, (0.12) 

which can be directly obtained from (0.5), (0.6), (0.7). The equality (0-12) and its generalization 
are of fundamental importance in the elaborated theory. 

The equality of the form AX — XB = C (where A, C arc square matrices ) is known as 
the Sylvester-Lyapunov equality (or as the Sylvester- Lyapunov equation, if it is considered as an 
equation with respect to X). 



^As usual, for the matrix M, denotes the transpose one. 
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Multiplying the right hand sides of the representations (0.9) and (0.11) term by term, we 
obtain the equality 

r{x)r-\y) = 1 + {x - y) eixl - A-p)-^S-\yI - A^^^e* (0.13) 

after some calculations. These calculations are based on the Sylvester-Lyapunov equality (0.12). 
The representations (0.9) and (0.11) are exactly what we need. They are said to be the system 
representations of the function r and r~-'-respectively. The representation (0.13) is said to be the 
joint system representation of the pair of (mutually inverses) functions r and r~^. 

We obtained the joint system realization (0.13) from the system representations (0.9) and 
(0.11). In its turn, the representations (0.9) and (0.11) are contained in (0.13): (0.9) is (0.13) for 
y = oo, (0.11) is (0.13) for x = oo. 

Now we derive the joint system representation (0.13) in a different way. Let k be the function 

ki.,y^= <^y-'(y)-' (0.14) 

x-y 

of two variables. Fixing y ^ ('P(r) U M{r) U oo), we consider k as a function of the variable x. 
This function is rational (with respect to x), with simple poles located at the points Xp G V{r), 
with the residues PXj,(y)' 

It is clear, that k{oo,y) = and that k{x,y) is holomorphic for x = y. Thus, k{x,y) admits the 
simple fraction expansion 

k{x,y)= Yl (0.16) 

l<p<n ^ 

In its turn, the residue PXp(y), considered as a function of y, is a rational one. Since r~^{Xp) = 0, 
the point Xp is not a pole of the function pXp, and the points fiq G J\f(r~^) arc the only poles of 
the function px^- All these poles are simple, with the residues hp^g = respAp(y) \y=ixq- 

_ 1 1 1 

As pXp{oo) = 0, the function px^, admits the simple fraction expansion 

l<9<n ^ ^1 

Substituting (0.18) into (0.16), we come to the formula 

k{x,y)= V -^^-K^-^—. (0.19) 

where hp^^ are defined by (0.17). The last formula can be presented in matricial form. Let us 
organize the numbers hp^q into the matrix 

-f^ = ||^p,g||l<p,g<n. (0.20) 

Thus, 

fc(x, y) = e{xl - Ar)-^HiyI - A^y^e*. (0.21) 
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and 

r{x)r-\y) = 1 + e{xl - Ap)-^H{yI - Aj^y^e*. (0.22) 
The comparison of the formulas (0.13) and (0.22) suggests us that H = S~^. The equaUty 

HS = I (0.23) 

may be verified starting from the formula 

The matrix identity (0.23) is equivalent to the system of scalar identities 

^ hp,gSg,t = il<p,t<n, py^t) (0.25) 

l<q<n 

and 

Yl Vs</,P = 1 (l<P<n). (0.26) 

l<q<n 

Prom (0.24) (and of course, from (0.5) and (0.17)) it follows that^ 

since r~^(Afe) = for all € 'P{r){= J\f{r~^)). Thus, (0.25) is verified. Analogously, the equality 
(0.26) may be verified using the formula 

Thus, the equality (0.23) holds. Finally, the formula (0.22) can be presented in the form (0.13), 
with 5" of the form (0.5). 

REMARK 0.1. By the way, we obtained the following rule for the inversion of the matrix S 
of the form (0.5), where Ai, A2, . . . , A„; /ii, 1x2, ■ ■ ■ , are pairwise different complex numbers 
(this matrix is known as the Cauchy matrix). Starting from these numbers, we have to construct 
the rational function r of the form (0.1) (with some c : c 7^ 0, 00). The entries hp^q of the inverse 
matrix H = are of the form (0.17). The expression (0.17) for the inverse matrix may 
be represented in the matricial form: 

= D,SD2, where D, = diag [ ^^Jy^^^^ ) > ^2 = diag . (0.29) 

In particular, we established the invertibility of the Cauchy matrix S. Of course, the expression 
(0.17) for the entries of the inverse matrix differs from the expression obtained in ([BGRl]) in 
form only. In ([BGRl]) (Lemma 6.3 there) this expression was obtained by means of the Kramer 
matrix inversion rule, using explicit expression for the Cauchy determinant. (The determinant 
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Ap — fiq fXq — At Xp — At \ Ap — fiq Xf — fJ'q 



5 



of the matrix S as well as its cofactors arc Cauchy determinants). The square of the Cauchy 
determinant can be calculated from (0.29): 

{detSf = n (r-'YiX,) . n r'in,). (0.30) 

REMARK 0.2. This derivation of the joint system representation can be extended to some 
classes of meromorphic functions. 

Obtaining "system" representations of "scalar" rational functions, wc made no use of anything 
that can not be generalized to the matricial case. Generalizations of the representations (0.9) and 
(0.11) to rational matrix functions have been done already. (See [S3], [GKLR], [BGRl], [BGR2], 
[BGRa]). 

The term "system representation" is related to the system theory. In this theory to each 
linear stationary time invariant dynamical system its transfer function is corresponded, which is 
matrix valued (and even operator valued) one. If the state space of the system is finite dimen- 
sional, then this transfer function is rational. And if the system is a SISO system (Single Input, 
Single Output), then its transfer function is scalar (complex vahicd). This transfer function is 
expressed in terms of the input, output and state space operators of the system. Inversely, each 
rational matrix function may be represented as the transfer function of an appropriate linear 
time invariant system with finite dimensional state space. Such a representation ^ is said to be 
the system realization or the system representation of the given rational function. 

The relationship between rational functions and linear systems can be exploited in both 
directions. We may apply results from matrix functions theory for study of linear systems. And 
we may use the system representation as a tool for study of rational matrix functions or as a 
tool to specify them. 

It should be mentioned that it was M. Livsic who discovered the characteristic matrix function 
of linear operators and operator colligations (=opcrators nodes). He was the first to relate 
invariant subspaces of a linear operator and factors of its characteristic matrix function, [L2], 
[L3], [LP]. In the system theory language, the characteristic function of an operator node is 
just the transfer function of the appropriate linear time invariant system. This was shown by 
M.S. Livsic in [L8], [L9], [BrL]. He has also related the characteristic matrix function with the 
scattering theory ([L5], [L6], [L9], [BrL]). 

The main goal of this paper is to present some basic results on the system representation of 
rational matrix function in a self-contained form. For the clarity of presentation, we restrict our 
consideration to the simple but important case of a rational matrix function in general position. 
The presentation is oriented to a "traditional" analyst. No previous knowledge in realization 
theory of matrix functions or its ideology is assumed. We leave detailed historical remarks for 
Section 5. 



1 . RATIONAL MATRIX FUNCTIONS IN GENERAL POSITION 

DEFINITION 1.1. Let M(z) be a k x k matrix function holomorphic in a punctured neighbor- 
hood of some point a (i.e. the point a is an isolated singularity of the function M). The point a 
is said to be a simple pole of the matrix function M if 

M{z) = ^+H{z), (1.1) 



It looks like the representation (0.9) 
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where Ma is a constant matrix and the matrix function H is holomorphic at the point a. The 
matrix Ma is said to be the residue of the matrix function M at the point a. 

DEFINITION 1.2. Let a point a be a simple pole of a k x k matrix function M and let the 
residue Ma of the function M at the point a be a matrix of rank one. As a matrix of rank one, 
the k X k matrix Ma can be factorized in the form 

Ma = fa- 9a (1-2) 

where fa is k vector-columns (i.e. k x I matrix), ga is k vector-row (i.e. 1 x k matrix), fa 7^ 
0,ga 7^ 0. The vectors fa and ga are said to be respectively the left semi-residual vector and the 
right semi-residual vector at the point a. The vectors fa and ga are defined uniquely, up to a 
constant factor: we can represent the number 1 in the form 1 = d- d"^ and then redistribute the 
factors d and d~^: 

fa^ fa- d; ga d~^ ■ ga. (1.3) 

We emphasize that the notions of left and right semi-residual vectors are defined only for a simple 
pole with residue of rank one. 

DEFINITION 1.3. A rational matrix function R (R € TZ{DJlk) ) is said to be a rational matrix 
function in general position if: 

1. The pole set V{R) and the zero set J\f{R) do not intersect; 

r{R)nJ\f{R) = 0; (1.4) 

2. All poles of the function R are simple, and the residues at these poles are matrices of rank 

one; 

3. All poles of the function R~^ are simple, and the residues at these poles are matrices of 
rank one; 

4. Both functions R and R~^ are holomorphic at the point z = 00. 



Let Rhe a rational matrix function in general position, with the pole set V{R) and the zero 
set J\f{R). For XeV{R), J^{R), let 

Rx = fx-gx, (/a7^0,5a7^0); R^ = f^,-g^, (/^^0,5m7^0) (1.5) 

be factorizations of the residue Rx of the function R at the point A and the residue R^ of the 
function R~^ at the point fx respectively, where /a,5a; /jU)5/i the appropriate semi-residual 
vectors. The additive expansions 

R(z) = Ri^)+ (1-6) 
xeP{R) 

R-\z) = R-\oo) + Yl ^ (1-7) 
ij.eAf{R) ^ 



can be rewritten in the form 



R{z)=R{oo)+ J2 /a-J^a'^"^' ^^-^^ 
xer{R) 
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R-\z) = R-Hoo) + J2 (1-9) 

LEMMA 1.1. I. Let a matrix function R be analytic in a punctured neighborhood of a point 
A € C, A be a simple pole of the function R and a holomorphicity point of the function R~^, 
and let Rx be the residue of the function R at the point A. Then the point A is a simple pole 
for the ^'logarithmic derivative" R' ■ R~^, and for the residue Px of this "logarithmic derivative" 
at this point the conditions 

Pi = -Px, rankPx = ranki?A (1.10) 

hold; in particular, 

tvacePx = — ranki?^- (l-H) 

11. Let a matrix function R~^ be analytic in a punctured neighborhood of a point fi £ C, fi 
be a simple pole of the function R~^ and a holomorphicity point of the function R, and let R^ 
be the residue of the function R^^ at the point Then the point jj, is a simple pole for the 

"logarithmic derivative" R' ■ R~^, and for the residue P^ of this "logarithmic derivative" at this 
point the conditions 

Pi = P^, rankP^ = ranki?^ (1.12) 

hold; in particular, 

traceP^ = ranki?^. (1-13) 



PROOF. We prove only statement I of Lemma. Statement II can be proved analogously. Let 

R{z) = + Ao + Ai{z-\) + A2{z - A)2 + . . . (1.14) 
z — A 

and 

R-\z) = Bo + Bi{z-X)+ B2{z -\f + - (1.15) 
be the Laurent expansions of the functions R and R~^ respectively. Then 

R!{z) = - . \.. +Ai + 2A2{z-\) + ... (1.16) 
{z - Aj 

Multiplying the Laurent expansions term by term, we obtain from (1.15) and (1.16) 

R'iz) ■ R-\z) = -j^^ - ^ + i-RB2 + A,Bo) + 0{z - A). (1.17) 

yZ — Aj Z — A 

Substituting the Laurent expansions (1.14), (1.15) into the identity R~^{z)R{z) = / and multi- 
plying these expansions term by term, we obtain 

BoRx = (1.18) 

B^Rx + BqAq = I. (1.19) 
Analogously, from the identity R{z)R~^{z) = I we derive 

RxBo = (1.20) 

RxBi + ^0^0 = /• (1-21) 



8 



(z-X) 

the point A is a simple pole for the function R'{z)R ^(z), with the residue 



Let us examine the expansion (1.17). According to (1.20), the term ^ vanishes. Thus, 



Pa = -i^A^i. (1.22) 
iProm (1.21) and (1.22) it follows that 7 + = AqBq. Hence, 

{I + Px)Px = (AoBo) ■ i-RxB,) = -Ao{BoRx)Bi. 

According to (1.18), BqRx = 0. Thus (/ + Pa)Pa = 0, i.e. P| = -Pa. Finally, 
PaPa = -{I- AoBo)Rx = -Pa + Ao{BoRx). Since BqRx = 0, 

Pa = -PaPa (1-23) 

^Prom (1.22) and (1.23) it follows that rankP^ = rankPA. □ 

REMARK 1.1. iProm (1.19) and (1.20) it follows that PaPiPa = Rx- Since Pi = (P~^)'(A) 
(see 1.15), it can be written as 

Rx{R-^y{X)Rx = Rx (AGP). (1.24) 

Analogously, we derive 

Rf,R'{li)Rf, = R^ ifieM). (1.25) 
We shall use the relations (1.24) and (1.25) in the following section. g 
LEMMA 1.2. For a rational matrix function R in general position, 

#P(P) = #AA(P) ( = n(P)). (1.26) 

PROOF. To prove the statement of Lemma, we have to consider, in some way or another, the 
"logarithmic derivative" R'{z) ■ R~^{z), (or R~^{z) • R'{z)). For a rational matrix function P in 
general position, its logarithmic derivative R' ■ R~^ admits the expansion 

P'(.).p-(.)= E E (1-27) 

Indeed, the (rational) function R'{z) ■ R^^{z) may have singularities only at the points of the 
set V{R) Uj\f{R). According to item 4 of Definition 1.3, the logarithmic derivative vanishes at 
the point oo; moreover, 

R'{z)-R-\z)=0{\z\-'^) {z^oo) (1.28) 

According to Lemma 1.1, all singularities of the logarithmic derivative are simple poles with 
residues Px and P^ of rank one (A € P(P), ^ G A/'(P) respectively). Thus, the expansion (1.27) 
holds, with tracePx = —1, traceP^ = 1. From (1.27) and (1.28) it follows that 

E Px+ E ^M = 0. (1.29) 

Because, according to Lemma 1.1, tracePA = —1, (A G P(P)), traceP^ = 1, (/x G M{R)), from 
(1.29) it follows that ^ (~1) + E 1 = 0- "^^^^ equality coincides with (1.26). Of course, 

xeP{R) ^l&^f{R) 

this result could be obtained using the operator version of Rouche theorem from [GS]. □ 
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DEFINITION 1.4. For a rational kxk matrix function R in general position, let us order its poles 
and ''zeros" somehow: P = {Ai, . . . , A„}, A/" = {/xi, . . . , /Xn}. (We remind that # P(R)=# N(R).) 
We introduce the n x n diagonal matrices 



A-p = diag(Ai, . . . , A„), Aj<^ = diag(/[Xi, . . . , fj,n), 

k X n matrices F-p, Fj^/ and nxk matrices G-p, G/^f: 

Ft = [f\i ) • • • ) /a„] , Fx = [/fti , • • • , //i„] , 



G-p 



9Xi 



9Xn 



Sfln 



(1.30) 

(1.31) 
(1.32) 



where /a , gx (Xj G V) are the left and right semi-residual vectors at the pole Xj of R, 
ffij^ dfij {f-i-j G A/") arc the left and right scmi-rcsidual vectors at the pole fij of R~^ . 

The matrices A-p and Aj^ are said to be the pole matrix and the zero matrix respectively for 
the matrix function R. 

The matrices Fp and Gp are said to be the left- and the right semi-residual matrices corre- 
sponding to the pole set V{R). 

The matrices Fj^ and Gj^f are said to be the left- and the right semi-residual matrices corre- 
sponding to the zero set J\f{R) . □ 

REMARK 1.2. It should be mentioned that if we order somehow the poles and the zeros, then 

the pole and the zero matrices Ap and Aj^ are defined uniquely, and the semi-residual matrices 
Fp, Gp,F_\f, G_\f are defined essentially uniquely, up to multiplication by diagonal matrices with 
non-zero diagonal entries: 



Fp — > Fp ■ Dp, Gp 



D 



where 



Gp, 



D-p = diag{dip. 



, d. 



'n,V> 



(d.p^O, j = l,2,...n). 



Dj^ = diag(di,Ar, • • ■ , dnX) (^iX 0, j = 1, 2, . . . n). 



(1.33) 
(1.34) 

(1.35) 
(1.36) 



This freedom in choice of these diagonal matrices Dp, Dj\f can be used to simplify some formulas. 
Of course, for given left semi-residual matrix F-p, the right semi-residual matrix G-p is determined 
uniquely; for given right semi-residual matrix Gp, the left semi-residual matrix Fp is determined 
uniquely, etc. □ 

It is clear that the additive expansions (1.8) and (1.9) may be rewritten in the matricial form 



R{z) = R{oc) + F-p ■ {zl - Ap)~^ ■ Gv, 
R-\z) = R-\oo)+FM-{zI-AMr^-GM. 



(1.37) 
(1.38) 



2. THE JOINT REPRESENTATION OF THE KERNELS ASSOCIATED 
WITH A RATIONAL MATRIX FUNCTION IN GENERAL POSITION. 

DEFINITION 2.1. Given a rational matrix function R of one variables, we associate with it 

r I 

two matrix function of two variables, K and K : 



R{x)R-\y)-I 



X 



y 



(2.1) 
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and 

^ X — y 

The function is said to be the right kernel associated with the function R. 
The function k'j^ is said to be the left kernel associated with the function R. g 

P(x) — P(y) 

REMARK 2.1. If P is a polynomial, the expression B(x,y) = is said to be the 

x-y 

Bezoutiant of the polynomial P. The expressions (2.1) and (2.2) look like a Bezoutiant. □ 

THEOREM 2.1. Let R be a rational k x k matrix function in general position; A-p and Aj^ be 
the pole and the zero matrices for R; F-p and Gj\f be the left pole and the right zero semiresidual 
matrices respectively. 
Then: 

1. The right kernel K^j^ is representable in the form 

K^x, y) = Frixl - A-py'H'-iyI - A^y^Gj^, (2.3) 

where is some n x n matrix. 

2. For given matrices ^ F-p and Gj\f, the matrix is defined uniquely: 



Tjr II II LT- ffAp f^lq ,^ .N 



3. The matrix = \\hl,J\ is invertible. 



The left version of this theorem holds as well. 



THEOREM 2.2. Let R be a rational k x k matrix function in general position; A-p and Aj^ be 
the pole and the zero matrices for R; Fj^ and G-p be the left zero and the right pole semiresidual 
matrices respectively. 
Then: 

1. The left kernel is representable in the form 

Ki{x, y) = Fm{xI - A^)-'H\yI - ApY^Gp, (2.5) 

where is some n x n matrix. 

2. For given matrices ^ Fj^ and Gp, the matrix = ||/ip,q||i<p,g<n is defined uniquely: 

= ll^p,gl|l<P,5<n, = ^I^I-Jh._ (2.6) 

3. The matrix is invertible. 

®We remind, that the matrices F-p, G-p, Fm, Gj^r are defined only up to transformations (1.33), (1-34) with 
arbitrary diagonal invertible matrices D-p, Dm- 
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DEFINITION 2.2. The matrices H'' and which appear in the representations (2.3) and 
(2.5) of the kernels and i^T^, are said to be the right core matrix and the left core matrix 
respectively. □ 

COROLLARY 2.1. Let R be a rational k x k matrix function in general position; A-p and 
Aj\f be the pole and the zero matrices for R; F-p, Fj^, G-p Gj^ are the appropriate semi-residual 

matrices; H^,H^ are the appropriate core matrices. 

Then the matrices R{x) ■ R~^{y) and R~^{x) ■ R{y) admit the representations 



R{x){R{y)r' = I+{x-y)F-p{xI-A-p)-'H'{yI-Aj^)-'Gj^, (2.7) 

{R{x))-^R{y) = I+{x-y)FM{xI-A^)-^H\yI-AvT^Gv. (2.8) 

Under the normalizing condition 

R{oc) = I, (2.9) 
the matrix functions R, R~^ themselves admit the representations 

R{z) = I - F'p{zl - A-pY^H^Gj^, (2.10) 

R-^{z) = I + FpH''izI - A^y^G^, (2.11) 

R{z) = I + FmH\zI-Ap)-^Gp, (2.12) 

R-\z) = I-Fm{zI-A^)-^H^Gp. (2.13) 



PROOF. The representations (2.7) and (2.8) are nothing more than the representations (2.3) 
and (2.5) rewritten in the terms of the functions R{x) ■ R~^{y) and R~^{x) ■ R{y). 
Letting y tend to oo in (2.7), we obtain (2.10); letting x tend to oo in (2.7), we obtain (2.11); 
Letting y tend to oo in (2.8), we obtain (2.12); letting x tend to oo in (2.8), we obtain (2.13). □ 

THEOREM 2.3. Let R be a rational kxk matrix function in general position; Ap and Aj^ be the 
pole and the zero matrices for R; Fp, Fj^, Gp G_\f are the appropriate semi-residual matrices; 
H^, are the appropriate core matrices. 

Then these matrices axe involved in the Sylvester-Lyapunov equalities: 

ApH^ - H'-Am = GpF^, (2.14) 
Aj^H^ - H^Ap = G^Fp. (2.15) 

PROOF. The matrices Gp>Fj^ and Gj^F-p are of the form 

Gp ■ F^f = \\gXj, ■ ffMj\l<p,q<n, Gj^Fp = \\g^p ■ fXg\\l<p,q<n- (2-16) 

The assertion of Theorem 2.3 follows from the explicit expressions (2.16), (2.4), (2.6), (1.30). g 

THEOREM 2.4. Let Rbe a rational kxk matrix function in general position, which satisfies the 
normalizing condition (2.9); Fp>, F^, Gp> Gj^ are the appropriate semi-residual matrices; H^, 
are the appropriate core matrices. 

Then these matrices axe involved in the equalities 

a). H'-Gm = -G-p; h). H^Gp = -Gm; c). FpH"- = F^; d). F^H' = Fp (2.17) 
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We give two proofs of Theorem 2.4. 

PROOF I. We compare the formulas (1.37) and (2.10). In the additive representation (1.37), 
for given left semi-residual matrix F-p, the right one G-p is determined uniquely. (See Remark 
1.2). Therefore, (2.17.a) holds. Analogously, comparing the formulas (1.37) and (2.12), we obtain 
(2.17.d). Comparing (1.38) and (2.11), (2.13), we obtain (2.17.c) and (2.17.b). □ 

PROOF II. We proof only the equality (2.17.a). The equality (2.17.c) can be proved analogously. 
(2.17.b) follows from (2.17.a) and (2.31), etc. The matricial equality (2.17.a) is the same that 
the system of n scalar equalities 

K,i9^. = -9X„ p=l, 2...,n. (2.18) 

l<q<n 

Since fx^ ^ 0, the last equality is equivalent to the equality 

f^p'i H ■ = -hv ■ 9Xp- (2.19) 

l<q<n 

Substituting (2.4) into (2.19), we come to equality 

• ( E ■ 9,. = -K ■ 9X., (2.20) 



l<q<n ~ ^1 



which is the same that the equality 



/Ap5A, • ( T^) = -/Ap5A,- (2.21) 



l<q<n P 

Since fx^gx^ = Rxp, Ugg^g = Rug, the last equality takes the form 



l<q<n 



(2.22) 



According to (1.7), 

Hence, (2.22) takes the form Rx^ ■ R~^{\p) = 0. According to (1.20), this equality is true. □ 



The proofs of theorems 2.1 and 2.2 are analogous. We will prove only the first one of them. 

PROOF of Theorem 2.1 . I. First of all, we obtain the representation (2.3). The main idea of 
the proof is to expand the function Kj^ of two variables into a double simple fraction series and 
then to interpret this expansion as the matricial equality (2.3). Actually, we derive not a double 
expansion of a function of two variables, but an iterated one. Let us fix a point y € C\M{R). For 
this fixed value y, we consider the kernel K'^[x, y) as a function of the variable x. This function 
is rational with respect to x, and ifj(x,j/) — > by x — > oo. It may have singularities only at 
the points of the set V{R) and at the point y where denominator x — y vanishes. Actually this 
function is holomorphic at the point x = y because the numerator vanishes at the point x = y 
as well. At each of the points Ap G 'P(-R) the function K'^{x,y) (considered as a function of x) 
either is holomorphic or has a simple pole with residue K^^{y) of rank one: 

Rx ■R'Hy) 

= \ (1 < P < n). (2.23) 
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Expanding the function K^{x,y) into the simple fraction sum, we obtain 

l<p<n ^ 

In its turn the residue K^^{y), considered as a function of y, is a rational function. It van- 
ishes at the point oo. This function is also holomorphic at the point Ap because the numerator 
Rxj, • R^^{y) vanishes at the point A^: the equality Rx^ ■ R~^{Xp) = is the same as the equality 
R\Bo = in (1.20). Thus, the only possible singularities of the function K^^(y) are the points 
(g = 1, . . . ,n) of the set J\f{R). These singularities are simple poles, with the residues Kp ^: 

K;,, = ^t^' (l<P,9<n)- (2-25) 
Expanding the function K^^ into the simple fraction sum, we obtain 

l<q<n y ^« 

Combining (2.24) and (2.26) (and transforming the iterated sum into the double sum), we obtain 
the double expansion 

l<g<n 

Substituting into (2.25) expressions (1.5), we get 

Kp,q = f>^p ■ K,q ' •S'm,' (2-28) 
where hp ^ are defined by (2.4). Thus, the expansion (2.27) takes the form 

l<p<n ^ > y f^q 

l<q<n 

where and g^^ are the left semi-residual vector at the pole Ap and the right semi-residual 
vector at the "zero" /iq respectively. 

The representation (2.3) is simply the representation (2.29) in the matricial form. The state- 
ment 1 of Theorem 2.1 is proved. 

Prom (2.27) it follows, that 

K;^^ = hm {x - Xp){y - fiq) ■ Kl{x,y). (2.30) 

Thus, the values Kp^q are determined from the kernel K'^{x,y) uniquely. From (2.28) it follows, 
that (for given fx^, g^^) the values hp^q are determined uniquely. The statement 2 of Theorem 
2.1 is proved. 

The statement 3 of Theorem 2.1 follows immediately from Theorem 2.5 below, where we not 
only prove the invertibility of matrices Hr and Hr, but also find their inverse matrices. □ 



THEOREM 2.5. Let R be a rational matrix function in general position, and the matrices 
H' , arc dcEned from it according to (2.4), (2.6), where {Ai, . . . Xn} = V{R), {/ni, . . . jin] = 
N{R)}, fx, gx, fn, gij, are appropriate semi-residual vectors. Then the equality holds: 

H'' = 1, ■H'' = 1 (2.31) 
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REMARK 2.2. If wc already know from somewhere that the matrix H is invcrtiblc, then we 
can easily deduce that {H^)~^ = . Indeed, multiplying the equation (2.15) by the matrices 
{H^) ^ from the both sides, from the right and from the left, and taking into account that 

(these are equalities (2.17.b) and (2.17.d)) , we obtain the equality 

A-piHY^ - {H^r^A^ = GvFm. (2.32) 

Thus, each of the matrices {H^)~^ and is the solution of the same Sylvester-Lyapunov equa- 
tion A-pX - XAj^ = G-pFj^. The condition V{R) riJ\f{R) = means that gaj, n cr^ ^ = 0- Under 
this condition, the solution X of the Sylvester-Lyapunov equation A-pX — XA^ = G-pFj^ is 
unique. Hence, {H^)~^ = ■ □ 



We give two proofs of Theorem 2.5. 

PROOF I. We prove only the first equality in (2.31). Multiplying (2.14) by the matrix from 
the right and (2.15) by the matrix from the left, we came to the equalities 



A-pH^'H^ - H'^AmH^ = GvFj^H^ 



and 



H'-Aj^H^ - H'-H^A-p = H^'Gj^F-p. 
Taking into account the equalities (2.17.a) and (2.17.d), we obtain 



and 



A'pH^'H^ - H^'Aj^H^ = G-pF-p 



H^AmH^ - H^'H^Ap = -GpFp. 



Adding two last equalities, we see that the matrices H^H^ and Ap commute: 

{H'' H^) Ap = Ap H"" hK 

Hence, 

{H'-H^) ip{Ap) = ip{Ap) {H'-H^). (2.33) 

for every function ip which is holomorphic on the spectrum cr^p of the matrix Ap. ^Prom (2.17.a) 
and (2.17.b) it follows that 

{H^ H^^ Gp = Gp. 

Multiplying this equality by (p{Ap) from the left and taking into account the commutational 
relation (2.33), we obtain that 



{H^'H^ - I) ip{Ap) Gp = 0. 



(2.34) 



Let us fix an index g G [1, . . . , n] and specify the function ip : (p{Xp) = Sp^g, p = 1, 2, . . . , n. 
By such choice of ip, tp{Ap) = diag[6i^q, S2,q, Sn,q]- Thus, 



ip{Ap) Gp 
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Therefore we obtain that rup^q g^^ =0 for 1 < p,q < n, where M =^ — /, M = 

\\f^P,q\h<p,<i<n- Since mp^q G C and g^^ is a non-zero vector row, rup^q = for sl\ 1 < p,q < n, 
i.e. M = 0. Hence, H^'H^ -1 = 0. □ 

PROOF II. We prove only the first equahty in (2.31). This matrix equahty is equivalent to the 
system of scalar equations 

E K,<i ■ K,P = 1' P = 1, 2, . . . , n. (2.35) 

l<g<n 

and 

XI K,q ■ = 0, p = 1, 2, . . . , n; s = 1, 2, . . . , n; p 7^ s. (2.36) 

l<9<n 

According to (2.4), (2.6) equality (2.35) means that 

(Ap-/x„)2 " ^ 



Because ^ 0, 7^ 0) the last equality is equivalent ^ to the equality 

(Ap - ^lq) 



/Ap • 2^ _ „ N2 = -JXp ■ gXp- (2.38) 



l<9<n 

or, what is the same^, to the equality 

fxpgxp ■ X] (x^'^^\2 ' f\9Xp = -fxpgxp- (2.39) 

l<q<n ^ ^ "^-^ 

Taking into account the factorization (1.5), we see, that the equality (2.39) is equivalent to the 
equality 

• E (X \ )2 ■ = -^A,. (2.40) 

In view of (1.7), 

Thus, the equality (2.40) takes the form 

Rx^-iR-'yiXp)-R^^ = Rx^. (2.42) 

According to Remark 1.1, the equality (2.42) holds. (See (1.24)). Thus, the equalities (2.35) are 
established. 

The equalities (2.36) can be established in the same way. According to (2.4) and (2.6), the 
equality (2.36) means that 

y i9XpUq)-igyxJ ^ _ 



^ If / is non-zero vector-columns, g is non-zero vector-row, then the equahty ci = C2, where ci, C2 are complex 
numbers, is equivalent to the equality fcig = fc2g. 



Here we use the associativity of the matrix multiplication. 
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Because /aj, 7^ 0, ^ 0, the last equality is equivalent (see the footnote ^) to the equality 

or (see the footnote , what is the same, to the equality 

h.9x, ■ E - Y^] ■ hs9x. =0, p ^ (2.45) 

Taking into account the factorization (1.5), we see, that the equality (2.45) is equivalent to the 
equality 



l<q<n 



In view of (1.7), 



E - r^) = - R~"i>^s), P + s. (2.47) 

Thus, the equality (2.46) takes the form 

Rx, ■ [r-\Xp) - R-\\s)) ■ Rx. = 0. (2.48) 

In view of (1.20) and (1.18), 

i?A^-i?~'(Ap) = 0, R-\Xs) ■ Rx, = 0. (2.49) 

Thus, (2.46) holds. The equalities (2.36) are established. □ 

The representations (2.7) and (2.8) are almost what we need. However, there is an essen- 
tial disadvantage in these representations: each one of them contains explicitly all the four 
semi-residual matrices. For example, the representation (2.7) contains explicitly not only the 
semi-residual matrices F-p and Gj\f (this is evident), but also the matrices Fx and G-p (see the 
expression (2.4) for the the right core matrix H'^). But the four semi-residual matrices (together 
with zero and pole locations) are over-determined data: the matrix function R is completely 
determined by two of those semi-residual matrices only. For example, from the additive repre- 
sentation (1.37) it follows, that (under the normalization R{oo) = I) the zero and pole matrices 
A-p, Ax together with the two semi-residual matrices F-p, G-p determine completely the ma- 
trix function R, and hence the other two semi-residual matrices Fx, Gx- Because of this, the 
semi-residual matrices Fx, Gp are (at least in principle) expressible in terms of the matrices 
Ap, Ax, Fp, Gx- Hence, we can hope to express the right core matrix in terms of the ma- 
trices Ap, Ax, Fp, Gx- Indeed, this can be done easily and explicitly: on the one hand, the left 
core matrix i?' is expressible in terms of the entries of the matrices Ap, Ax, Fp, Gx only (see 
(2.6); on the other hand, H'' = {H^y^ (see (2.31)). 

This suggests us that it may be reasonable to use the inverse matrices and 5' instead of 
the matrices^ and H^: 



*We remind, that, according to Theorem 2.3, the matrices and are mutually inverse, and hence, 



invertible. 
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The equalities (2.17) and (2.14), (2.15) can be rewritten in term of the matrices S"" and S . 
Namely, equalities (2.17) take the form 

a).Gj^ = -S''Gr; h). Gv = S^Gm; c). F-p = FmS'; d). Fm = F-pSK (2.51) 

Multiplying the equality (2.14) by the matrix (H''')^^ from the right and from the left and 
taking into account equalities (2. 17. a) and (2.17.c), we transform (2.14) to the form 

Aj^S'' - S^'A'p = Gj^F-p. (2.52) 
Analogously, from (2.15) and (2.17.b), (2.17.d) we derive the equality 

ArS^ - S^Am = GvFm- (2.53) 
Thus, the matrices S'^ and are solutions of the Sylvester-Lyapunov equations 

Aj^X - XA-p = Gj^F-p. (2.54) 

ApX - XAj^ = GpF^. (2.55) 

respectively. 

Now we change our point of view and define the matrices and S*' as solutions of Sylvester- 
Lyapunov equations (but not as the matrices inverse to the core matrices and H^; see (2.50)). 

DEFINITION 2.3. Let i? be a rational matrix function in general position, Ap and Aj\f be its 
pole and zero matrices, Fp, Gp>, Fj\f, Gj\f be its appropriate semi-residual matrices. 

1. The matrices and which arc the solutions of the Sylvester-Lyapunov equations (2.54) 
and (2.55), are said to be the right zero-pole coupling matrix and the left zero-pole coupling 
matrix respectively. 

2. The relations (2.51) are said to be the zero-pole coupling relations. 



REMARK 2.3. Since the spectra of the matrices Ap and Aj\f do not intersect, the Sylvester- 
Lyapunov equations (2.55) and (2.54) are uniquely solvable. (However, as the matrices Ap 
and Aj^ are diagonal, the solvability of these equations as well as the uniqueness is obvious). 
Moreover, it is possible to obtain the explicit expressions for the matrices and S'^ from (2.55) 
and (2.54): 

r). S"" = \\Sp^g\\l<p,q<n, = ^ , 1). 5' = \\Sp^q\\l<p,q<n, = ^ ^ . (2.56) 

H'p ^p h^q 

(Actually, we derived the Sylvester-Lyapunov equations from the explicit expressions for the 
matrices which we interpret now as solutions of these equations). □ 

According to (2.50), the relations (2.31) can be rewritten in the form 

S'' ■ = ■ S"" = I. (2.57) 

We may also refer directly to the equalities (2.37) and (2.43): these equalities mean that the 
matrices (2.56.r) and (2.56.1) are mutually inverse. 
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REMARK 2.4. The representations (2.3) and (2.5) may be rewritten in terms of the matrices 
and (instead of the matrices and ii"'): 

(i?(x)(i?(y))-i = I + ix-y)Fr{xI-Ar)-' ■{sn-^■{yI-A^ry^G^f, (2.58) 
{R{x)r'R{y) = I + {x-y)FM{xI-A^)-^-{S')-^-{yI-A-p)-^G-p. (2.59) 

Under the normahzing condition (2.9), the matrix functions R,R~^ themselves admits the rep- 
resentations 



R{z) = I - F-p{zI - A-pJ-^S'-y^GM, (2.60) 

R-\z) = I + F-p{S'')-\zI-Am)-^Gm, (2.61) 

R{z) = I + F^/{S^y\zI - Av)-^Gv, (2.62) 

R-\z) = I-FM{zI-AMy\S^y^Gv. (2.63) 



These formulas may be obtained of from (2.58), (2.59), letting x or y tend to oo there, or from 
(2.10) - (2.13), rewriting them in terms of the matrices S"',S'' (instead of the matrices H^,H^). 

□ 

REMARK 2.5. Of course, we may obtain the zero-pole coupling relations (2.51) comparing the 
representations (1.37), (1.38) and (2.10) - (2.13). □ 

REMARK 2.6. The semi-residual vectors are defined not completely uniquely, but up to 
transformations (1.33) and (1.34) only. If the semi-residual vectors are transformed according 
to (1.33) and (1.34), the right hand sides of the Sylvester-Lyapunov equations (2.54) and (2.55) 
are transformed as: 

GmFv ^ {Dj^y^ ■ GmF-p ■ Dj,; G-pFm ^ {D-py^ ■ GvFm ■ Dm (2.64) 

The solutions S'^ and of the Sylvester-Lyapunov equations (2.54) and (2.55) are transformed 
as: 

s"- ^ {D^y^ ■ s"- ■ Dv- ^ {Dpy^ ■ ■ d^. (2.65) 

Of course, the expressions (2.58), (2.59) (for R{x){R(y))~^ and (i?(x))~^i?(y)) are invariant 
with respect to the transformations (1.33), (1.34), (2.65) (of the semi-residual and the zero-pole 
coupling matrices). □ 

The representations (2.58) and (2.59) are exactly what wc need. Obtaining them is one of 
the main goal of this paper. Therefore wc choose a special name for this representation: 

DEFINITION 2.4. Let R be a rational matrix function in general position, A-p, Aj\/ be its 
pole and zero matrices, F-p , G-p , Fj^/ , Gx be appropriate semi-residual matrices, S^, iS' be the 
solutions of the Sylvester-Lyapunov equations (2.54) and (2.55) respectively The formulas (2.58) 
and (2.59) are said to be the right joint system representation of the pair R, R~^ and the left joint 
system representation of the pair R, R^^ respectively. 

The formulas (2.60), (2.61), (2.62), (2.63) (which can be obtained from (2.58), (2.59) by 
passage to the limit) are said to be the right system representation of the function R, the right 
system representation of the function i?^, the left system representation of the function R, the left 
system representation of the function R~^ respectively. 

REMARK 2.7. The terminology is motivated by the so-called system theory or, in more detail, 
by the theory of linear time invariant dynamical system. In this theory, all the objects such 
as the zero and pole matrices, the semi-residual matrices, the zero-pole coupling matrices are 
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interpreted from the point of view of dynamical systems. This interpretation does not play any 
role in our considerations. Wc need the joint system representations as a tool to introduce a 
convenient coordinates in the set of all rational matrix functions (in general position). g 

REMARK 2.8. In the realization theory one obtains formulas like (2.60) - (2.63) for matrix 
functions R and considered individually. In the representations (2.58), (2.59) the matrix 
function R, R~^ are considered jointly. This is the reason for using the terminology joint system 
representation. g 

Now we summarize the results of this section and formulate 

THEOREM 2.6. Let R be a rational matrix function in general position, A-p and Aj\f be its 
pole and zero matrices, Fp, G-p be its left and right pole semi-residual matrices, Fj\f, Gj\f be its 
left and right zero semi-residual matrices. Then: 

1. The matrices S"", which is a solution of the Sylvester-Lyapunov equation Aj^j-X — XA-p = 
Gj\fFp, and S'', which is a solution of the Sylvester-Lyapunov equation A-pX — XAj^ = 
GvFm: 

^\ qr_\\„r II r _ 9tip ' f\q ,n nl _ \i I || J _ 9\p ' f^lq (c) r^c) 

I). J — ||''p,g||l<p,g<n) ~ ,, _ \ ' '■)■ i~> — \\^p,q\\l<p,q<n: '^p^q — _ j {^-OV) 

are mutually inverse, i.e. the equahties 

S'' -S^ = 1, 5' • 5^ = / (2.57) 

hold. In particular, the matrices S^' and S"' are invertible. 

2. The matrix function R{x){R{y))~^ can be recovered from the data Ap, Aj^, Fp, Gj^ by 
the formula (2.58) 

iRix){R{y))-' =I+{x-y)Fp (xl - A-p)-' ■ {S')-' ■ {yl - A^y^G^, (2.58) 
i.e. the right joint system representation holds. 

3. The matrix function {R{x))~^ R{y) can be recovered from the data Ap, Aj\f, Fj\f, Gp by 
the formula (2.59), 

{Rix))-'Riy) = I+(x-y)FM (xI - A^y' ■ (S'y' ■ {yl - Apy^Gp. (2.59) 
i.e. the left joint system representation holds. 

4. Under the normalizing condition (2.9), the matrix functions R, R~^ themselves admit the 
representations 

R{z) = I-Fp{zI- Ary\S'y^Gu, (2.60) 
R-\z) =I + Fp [S'-y^zI - Aj^y^Gj^, (2.61) 
R{z) = 1 + Fj^{S^y\zI - A-py^Gv; (2.62) 

R-\z) =1- Fm{zI - A^y^S^y^Gp. (2.63) 

5. The zero-pole coupling relations hold: 

a). Gm = -S^Gp; b). Gp = -S^Gj^; c). Fp = F^S""; d). F^ = FpSK (2.51) 

6. If the representations (2.7) and (2.8) hold with some matrices and , then 
jjr ^ (S'-y^, = (S^y^ of necessity 
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REMARK 2.9. The representation (2.58) allows us to recover the matrix function R{x){R{y))~^ 
from the left pole- and the right zero- semi-residual matrices F-p and Gj\f. whereas the representa- 
tion (2.59) allows us to recover the matrix function {R{x))~^ R{y) from the right pole- and the left 
zero- semi-residual matrices Fa/ and G-p. However, sometimes one have needs for some "hybrid " 
formulas which allow to recover the matrix function R{x){R{y))~^ from the right pole- and the 
left zero- semi-residual matrices F_\f and G-p and the matrix function {R{x))~^ R{y) from the left 
pole- and the right zero- semi-residual matrices F-p and Gj^/. Such formulas can be easily derived 
from the joint system realization formulas (2.58), (2.59) combined with the zero-pole coupling 
relations (2.51). These "hybrid" formulas are of the form: 

Rix) ■ (i?(y))-^ = I-{x-y)F^f (5')-^ (xl - A-p)'^ (yl - A^f)-^ {S'y'Gp, (2.66) 
{R{x))-' ■ R{y) = I-{x-y)F-p {S^-' {xl - A^Y^ {yl - A^Y^ {S^Y^G^. (2.67) 

The matrix S*' can be calculated from the data: Fj^^ Gp, A-p, Af^; the matrix can be calcu- 
lated from the data: F-p, Gn, A-p, A_\f. □ 

3. FROM THE CHAIN IDENTITY TO THE SYLVESTER - LYAPUNOV 
EQUATION AND BACK. 

The consideration of this item are concentrated around of the so-called chain identity. Let 
us give a number of definitions. 

Let T{. , .) he a k X k matrix function of two complex variables, with domain of definition ^ 
Vt,Vt eCxC, T:Vt^ Mk- 

DEFINITION 3.1. A function T{., .) of two variables is said to satisfy the chain identity if 

T{x,y)-T{y,z) = T{x,z) (3.1) 

for every x, y, z for which {x, y) G Vt and {y, z) G Vt- (In particular, (x, z) must belong to Vt, 
if {x, y) e Vt and {y, z) G Vt). 

DEFINITION 3.2. A function T{. , .) of two variables is said to satisfy the diagonal unity 
identity if 

T{x,x)=I (3.2) 
for every point x for which (x, x) belongs to the domain of definition Vt of the function T. 

A class of function T satisfying both the chain identity and the diagonal unity identity can 

be constructed in the following way. 

DEFINITION 3.3. Let $ and be a k x k matrix functions of one variables with domains 
of definition P$ and V^-i respectively, V^ G C, V^-i G C. Let us define the matrix function 
T$ of two variables by the equality 

T^ix,y)''^'<^{x)-<^-\y) (3.3) 

with domain of definition 

Vt = P$ X V^-i, i.e. (^{x, y) G Pt) ^ ({x G P$) & (y G P$-i)) . (3.4) 
The function r$ is said to be the chain function generated by the function 

LEMMA 3.1. Let $ and be a k x k matrix functions of one variables with domains of 
definition V^ and V^-i respectively, V^ G C, V^-i G C. Let r$ be the chain function generated 

^We recall that C = C U oo is the extended complex plane. 
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by the function If the matrix functions $ and ^~ are mutually inverse, i.e. the identities 
^{x) ■ <^>^^(x) = $~^(x) • ^{x) = I hold for aU a; € P$ n V^-i, then for the matrix function T 
both chain identity and diagonal unity identity hold. 

REMARK 3.1. Of course, Lemma 3.1 is reach in content only under condition 2?$ fi V^-i ^ 0. 
If this condition fails then the values T{x, x) and T{x, y) ■ T{y, z) are defined for the empty set 
of arguments. 

PROOF of Lemma 3.1. The diagonal unity identity expressed that the functions $ and are 
mutually inverse. The chain identity for the function T<j> is the consequence of two facts: 1). The 
function $ and ^>~^ are mutually inverse; 2). The matricial multiplication is associative. g 

It turns out that each function T of two variables satisfying both the chain identity and the 
diagonal unity identity is of the form r$ for some function $ of one variable. 

THEOREM 3.1. Let T be a k x k matrix function of two variables, which domain of definition 
Vt eCxC is of the form Vt = Vi x 1)2, where Pi G C and D2 € C, with Vir\V2^ 0. If for 
the function T both chain identity and diagonal unity identity are satisfied, then the function 
T is of the form T = T$, (i.e. T{x, y) = • ), where $ and are mutually inverse 

k X k matrix functions of one variable, with X>$ = Vi and T>^-i = X>2- 

PROOF. Let us fix an arbitrary point a belonging to the set X>i fl I>2- (We will call this point a 
the distinguished point.) Let us define now 

V^'^^'Vi, ^x)''^'T{x,ay, Vl^^^V2, = r(a,2/)- (3-5) 

The functions $ and arc mutually inverse: this follows from the chain and diagonal unity 
identities. The equality T(x,y) = ^>(x) • ^>^^(y) is the chain identity written down for the triple 
of the points x, a, y. In addition, we note that this function $ satisfies the normalizing condition 
$(a)=/. □ 

DEFINITION 3.4. Let R be a rational k x k matrix function of one variables, dctR ^ 0, and 
is the inverse (in the commonly accepted sense) matrix function; the domain of definition 
Vr is the set of holomorphicity of the function R; the domain of definition is the set of 
holomorphicity of the function (In other words, ^° = C \ V{R), = C \ N{R) ). We 

r I 

associate with the function R two functions of two variables, {x,y) and T^{x,y): 

a). T^ix,y) R(x) ■ R~\yy, b). r^(x,y) R-\x) ■ Riy). (3.6) 

r I 

The function is said to be the right chain function generated by R. The function T^ is said 
to be the left chain function generated by R. 

REMARK 3.2. It is clear that that the function T^ is the function of the form (in the sense 
of the Definition 3.3) for ^ = R, and the function T is the function r$ for $ = R~^. Thus, the 
right chain function generated by R is the left chain function generated by R~ : 

r;(x,y)=T^_i(x,y). (3.7) 



^°We recall that P{R) is the pole set of the function R, Af{R) is the zero set of the function R, i.e. the pole set 
of the function R~^. 
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^Prom Lemma 3.1 and from Remark 3.2 it follows 

COROLLARY 3.1. Let R be a matrix function, dctR ^ 0. Then botli matrix functions {x, y) 
and T^^{x,y) (see (3.6)) satisfy the chain identity and the diagonal unity identity: 

T; {x,y)- T; {y,z) = T; (x,zy, r; {x,x) = I. (3.8) 

(z, y) . tI (y, x) = (z, x) ; (x, x) ^ I. (3.9) 



^Prom (2.1) and (2.2) it is clear that the right and left chain functions and r^_i can be 
expressed from the right and left kernels: 

T;(x,y) = 7 + (x-y)-<(ar,y); (3.10) 
Tl(x,y) = I+(x~y)-Kl{x,y). (3.11) 

Substituting in the chain identity (3.8) the expression (3.10) for T^{x,y), we obtain an identity 



for the right kernel 

(/ + (x - y)Klix, y)) ■ (/ + (y - z)K'^iy, z)) = I + (x - z)Kl(x, z). (3.12) 
Removing the parentheses, we obtain the identity 

(x - y)Kl^(x,y) + {y - z)Kl{y,z) - (x - z)Kl{x,z) = -(x - y){y - z)K''^(x,y)Kl^(y,z). 
Dividing on (x — y){y — z), we come to the identity 

K^^(x, y) - K^ix, z) Kl{x, z) - Kl{y, z) 



y-z x-y 



= -K{x,y)K{y,z). (3.13) 



Assume now that the kernel admits the representation (2.3), with some k x n matrix F-p 
nxk matrix Gj\f and n x n matrices A-p,Aj\f,H^. Substituting the expressions (2.3) into (3.11), 
we come to the identity 

. - ^v)-' - - ^'■)-' .H-.i.I- Au)-' ■ Gu ^ (3.14) 

X-y 

- Fr ■ (xl - Ar)-' ■ ■ [yl - A^y' ■ GmFv ■ [yl - AvY^ ■ ■ {zl - A^Y^ ■ Gm- 
Using the Hilbert identity for resolvents, we come to the identity 

Fv ■ (xI - Ar)-^ ■ ■ (yl - Au)-^ ■ (zI - A^Y^ ■ Gm - 

F-p ■ (xl - A-p)-^ ■ (yl - Ap)-^ ■ H"- ■ (zl - Aj^)'^ ■ Gj^ = (3.15) 
Fp ■ (xl - ApY^ ■ ■ (yl - AmY^ ■ GmFp ■ (yl - ApY' " ■ (zl - A^fY^ ■ Gm- 

Putting the common factors outside the parentheses, we obtain 

Fp ■ (xl - ApY^ ■M-(zI- Aj^Y^ -Gn =0, (3.16) 
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where 

M = H'-{yI- AmT^ - {yl - At.)-^ -H'-H'- {yl - A^y^ ■ GuF-p ■ {yl - Av)'^ ■ (3.17) 



Assume moreover that i? is a rational matrix function in general position. Let A-p and Aj\f 
be its pole and zero matrices (i.e. these matrices are of the form (1.30), where all the numbers 
\i, . . . , Xn] Hi, ■ ■ ■ , IJ'n are pairwise different), and F-p and G_\f be the pole left semi-residual 
and the zero right semi-residual matrices (in particular, they have the form (1.31), (1.32) where 
fxp and g^^ are non-zero k vector-columns). 

^Prom (3.16) it follows that 

F-p ■ ^{A-p) ■ M ■ ^(Aaa) Gm =0, (3.18) 

where 99 and ^ are arbitrary functions which are analytic on the spectra of the matrices Ap and 
Aj^ respectively. 

Let us fix two indices, p G [1, . . . , n] and q . . . ,n\ and specify two functions ip and ^: 

^{^0) = ^j,P i^in) = U = i> 2, • , n). 

Then 

(p{Ap) = diag ((5i,p, . . . , 6n,p), ^{Af/) = diag ((5i,q, . . . , 5n,q). (3.19) 
For such functions (p and ip, (3.18) becomes the form: 

hp ■ mp,g ■ g^^ = 0, (3.20) 

where rUp^q is p, q entry of the matrix M. Since fx^ and g^^ are non-zero vector-column and 
vector-row, rup^q = 0. Thus, M = 0, or, 

■ {yl - Au)-^ - {yl - A-pY^ ■ - ■ {yl - Au)'^ ■ G^F-p ■ {yl - Ap)'' ■H' = Q. (3.21) 
The "left" version of the equation (3.21) has the form 

■ {yl - Ap)-^ - {yl - A^)-^ ■ - ■ {yl - Ap)-^ ■ G-pFj^ ■ {yl - Aj^)-^ ■H^ = 0. (3.22) 

To obtain (3.22), we have to use the chain identity (3.9), the expression (3.11) for the chain 
function Tjj in terms of its associated kernel Kj^, and the representation (2.5) of this kernel. 
However, we may just replace in (3.21): 

H^^H^; Fp^Fm; Gm^Gv- A^^A-p; A-p ^ Aj^. 

Let us examine the Laurent expansion (with respect to y) of the function in the left hand side 
of (3.21): 

■ {yl - Am)-^ - {yl - Ap)-^ -H'-H'- {yl - AmY^ ■ Gj^Fp ■ {yl - Ap)-' ■ = 

{H'Am - ApH' -H^Gm- FpH')y-^ + 0{y-^) {y ^ 00) 

In view of (3.21), the Sylvester-Lyapunov equality holds: 

H^'Aj^ - ApH'' = H^'Gm ■ FpH"". (3.23) 
Analogously, from (3.22) we can derive the equality 

H^A-p - Aj^H^ = H^Gp ■ FmHK (3.24) 
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According to Theorem 2.1 (actually, according to Theorem 2.5), the core matrix (see Defi- 
nition 2.2) is invertible. (The first proof of the Theorem 2.5 is based on the Sylvester-Lyapunov 
equalities (2.14), (2.15), but the second one is independent from them). Multiplying (3.21) by 
the matrix {H^)~^ from the right and by the matrix (iJ'")"^ from the left and denoting, as 

before (see (2.50)), 5^ = (ii"0~ > we come to the equality (2.52). ^Prom (3.22) we can derive 
the equality (2.53) in the same way as we already derived the equality (2.52) from (3.21). 

Thus, we obtained the equalities (2.52) and (2.53) in two different ways. The first one is 
based on the explicit expressions (2.4) and (2.6) for the core matrices and if. This method 
uses essentially the specific character of a rational matrix function in general position. The 
second method works for much more broad classes of rational matrix functions. Up to certain 
point, the method works for arbitrary rational matrix function. First of all, this method uses the 
chain identities (3.8) and (3.9). These identities are evidently true for arbitrary matrix functions 
R which are non-degenerate (i.e. deti? ^ 0). Then we use the representations (2.3), (2.5) for 
the kernels, associated with R. However, we nowhere use that the pole and zero matrices A-p 
and Aj\/ are diagonal, with disjoint simple spectra. Actually, we have obtained the equality 
(3.18) for amj rational matrix function R which associated kernel K'^^ admits the representation 
(2.1) with arbitrary A-p, Aj^, F-p,Gj\f and H' (or, what is the same, for any rational matrix 
function R, such that the function R{x) ■ {R{y))~^ admits the representation (2.7) with arbitrary 
A'p,Aj\/,F'p,Gj\/ and H'^). Then we have to conclude from (3.18) (under the assvmiption that 
this equality holds for arbitrary functions (p, tp which are holomorphic on the spectra of A-p and 
A-p respectively), that (3.21) holds. 

DEFINITION 3.5. (The row- version) : Let T be a k x n matrix (k rows, n columns), and A 
be n X n matrix. The pair (F, A) is said to be obstrollable^^ if the linear span of the set of 
n-vector-rows {vT{XI — A)~^}, where v runs over the space of all k-vcctor-rows and X runs 
over C\aA ( cta is the spectrum of A ), coincides with the whole space C" (of all n-vector-rows). 
The equivalent definitions: 

I. The pair (F, ^4) is said to be obstrollable, if the linear span of the set of vectors {v TA"^}, where 
V runs over of all k vector-rows and m runs over the set N of aiJ natural numbers, coincides 
with the whole space C" (of all n-vcctor-rows) . 

II. The pair (F,^) is said to be obstrollable, if the linear span of the set of vectors {vTip{A)}, 
where v runs over C'' of all k vector-rows and (p runs over the set of all functions holomorphic 
on GA, coincides with the whole space C" (of all n-vector-rows). □ 

DEFINITION 3.5 (The column- version): Let F be a n x k matrix (n rows, k columns), and b be 
n X n matrix, the pair (B, F) is said to be obstrollable, if the linear span of the set of n-vector- 
columns {{/jlI — B)~^rv} , where v runs over the space C*^ of aJi k-vector-columns and /x runs over 
C\aB ( CFB is the spectrum of B), coincides with the whole space C" (of all n-vector-columns). 
The equivalent definitions: 

I. The pair [B, F) is said to be obstrollable, if the linear span of the set of vectors {TB"^v}, 
where v runs over C'' of all k vector-columns and m runs over the set N of aiJ natural numbers, 
coincides with the whole space C" (of all n-vector-columns). 

II. Tlic pair {B, F) is said to be obstrollable, if the linear span of the set of vectors {f{B) Tv}, 
where v runs over of all k vector-columns and if runs over the set of all functions holomorphic 
on gb, coincides with the whole space (of all n-vector-columns). □ 

COMMENT TO TERMINOLOGY: This terminology is motivated by the system theory. (In 
^'^The word obstrollable is a mixture of the words observable and controllable. 
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more detail, by the theory of hnear time invariant dynamical systems). In this theory structures 
like F(XI — A)^^ and (fil — B)^^G appear, where F and G are k x n and n x k matrices and 
A, B are n x n matrix, and usually n is much bigger then A; (n » fe). If F is interpreted as the 
input operator of the system and G is interpreted as its output operator, then the notions of 
controllability of the pair {F.A) and the notion of observability of the pair {B, G) arc intro- 
duced. If G is interpreted as the input operator of the system and F is interpreted as its output 
operator, then the notions of controllability of the pair (S, G) and the notion of observability of 
the pair {F,A) are introduced. Structures analogous to the structure F{\I — A)~^ , {jil — B)~^G 
appear in the system representation problems as well. We need to formulate the property which 
is controllability if F is interpreted as the input operator, and is observability, if T is interpreted 
as the output operator. However, we would not like to give the preference to one of two: in- or 
out- interpretations of the matrix V. Because this, we choose the "neutral" term obstrollability. 

□ 

LEMMA 3.2. Let F,G he k x n and n x k matrices, and A, B, M are nxn matrices. Assume 
that 

1. F{XI - A)-^M{fiI - B)-^G = (VA € C \ (7^, V// e C \ as). 

2. The pairs {F,A) and {B,G) are obstroUable. 

Then M = 0. 

PROOF. The proof follows immediately from the definition of obstrollability. q 

Thus, the foregoing reasonings ( the reasoning of this section and the reasoning used for the 
first proof of Theorem 2.5 ) prove the following 

THEOREM 3.2. I). Let R be a rational function such tLat^^ the chain matrix function 
R{x){R{y))^^ admits the representation of the form (2.7) some k x n matrix F-p, n x k matrix 
Gj\f and nxn matrices A-p, Aj^ and . 

If the pairs {F-p, A-p) and {Aj^j-, G//) are obstroUable, then the the matrix satisfy the 
equahty (3.23); 

II) . Let R be a rational function such tiiat^^ the chain matrix function {R{x))~^ R{y) admits 
the representation of the form (2.8) some kxn matrix Fj\f, nxk matrix Gp and nxn matrices 
Ap,Aj^ and HK 

If the pairs {Fj^j-, Aj^j-) and {Ap, Gp) are obstroUable, then the the matrix satisfy the 
equality (3.24); 

III) . Let R be a rational function such that the chain matrix functions R{x){R{y))~^ and 
{R{x))^^R{y) admit the representations of the form (2.7) and (2.8) respectively, with some 
kxn matrices Fp, F^, sonic nxk matrices Gj\f, Gp and some nxn matrices Ap, Aj^, H^, HK 

If all four pairs {Fp, Ap), {Ax, G_\f), {F_\f, A_\f) and {Ap, Gp) are obstroUable, and if 
moreover the coupling relations (2.17) hold, then the matrices H^,H^ are mutually inverse: 
H"" ■ = I,H^ ■ H'' = I, and for their inverse matrices S'' = {H'')-^ and = (H^)'^ the 
equalities (2.52), (2.53) holds, i.e. the matrices and are solutions of the Sylvester-Lyapunov 
equations (2.54) and (2.55). 

THEOREM 3.3. Let R be a rational function such that the chain function R{x){R{y))~^ admits 

We assume neither the matrices A-p, Aj^r are diagonal, nor their spectra are simple, nor their spectra are 
non-intersecting. We also don't assume a'priori that the matrices H^,H' are invertible. 
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the representation of the form (2.58) with some k x n and n x k matrices F-p, Gj\f, n x n 
matrices A-p and Aj^ and some invertible nxn matrix S'^' . If the pairs {F-p, A-p) and {Aj\f, Gj\f) 
are obstrollable, then the equality (2.52) holds, i.e. the matrix S'' is a solution of the Sylvester- 
Lyapunov equation (2.54). 

The "left" version of this theorem holds as well. 

THEOREM 3.4. Let R be a rational function such that the chain function {R{x))~^R{y) 
admits the representation of the form (2.59) with^^ some k x n and n x k matrices F_\f, Gp, 
nxn matrices Ap and Aj^ and some invertible nxn matrix If the pairs {Fj^, ^Af) ^.nd 
{Ap, Gp>) are obstrollable, then the equality (2.53) holds, i.e. the matrix is a solution of the 
Sylvester-Lyapunov equation (2.55). 

LEMMA 3.3. Let F = (/i, . . . , /„) be a k x n matrix (i.e. f\, ... , fn are k vector-columns), 
and A = diag(ai, . . . , a„) be a diagonal matrix with simple spectrum (i.e. all diagonal entries 
ai, . . . a„ are pairwise different). If no column fi, ... , fn of the matrix F is equal to zero, then 
the pair (F, A ) is obstrollable. 

PROOF. Actually, the proof of the statement was already done (See how we obtained the 
equality (3.20)). □ 

This lemma (together with Theorem 2.1) shows that Theorem 3.2 and Theorem 3.3 are 
applicable to rational matrix functions in general position. 

Now wc "inverse" our reasonings leading from a chain identity to a Sylvester-Lyapunov 
equation. Let F, G he k x n and n x k matrices. A, B he n x n matrices, with spectra and 
aB, and S be an invertible nxn matrix for which the equality 

BS-SA = GF (3.25) 

holds. Evidently, this equality is equivalent to the identity (with respect to y G C): 

SiyI -A)- (yl - B)S = GF 

Multiplying the last identity by the matrix S~^{yl — B)~^ from the left and by the matrix 
{yl — A)~^S~^ from the right, we come to the identity 

S-^iyl - B)-^ - {yl - A)-^S-^ - S-\yI - B)-^G ■ F{yl - A^^S'^ = 0. 

(This is nothing more than the equality M = 0, where M is defined by (3.17)). Multiplying the 
last identity by the matrix F{xl — A)~^ from the left and by the matrix (zl — B)~^G from the 
right, we come to the identity (with respect toxSC, y € C, z G C): 

F{xl - A)'^S-^ ■ {yl - B)-^{zI - B)-^ ■ G- 
-F ■ {xl - A)-\yI - A)-^ ■ S-\zI - B)-^G- 
F{xl - A)-^S-\yI - B)-^Gx 
xF{yI - A)-^S-\zI - B)-^G =0. 

Using Hilbert identities 

x-y 

{yl - B)-^ - {zl - B)-^ 



{yl - B)-\zI - B)-^ = 



y- z 



^^We assume neither the matrices A-p, are diagonal, nor their spectra are simple, not their spectra are 
non-intersecting. Of course, the invertibility of the matrix S"" imposes implicitly some restrictions on the matrices 
A-p, Aj^, Fv, Gm- 
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we reduce the last identity to the form 

K{x,z)-K{jj,z) K{.r,y)-K{x,z) 



= K{x,y) ■ K{y,z) 



x-y y-z 
or, what is the same, to the identity 

{x - y) K{x, y) + {y ~ z) K{y, z) - {x - z) K{x, z) = -(x - y){y - z) K{x, z). (3.26) 

where now the function K{(,r)) of two variables is defined by the formula 

K{x, y) = F{xl - A)-^ ■ S-^ ■ {yl - B)-^G. (3.27) 

The identity (3.26) may be rewritten in the form a chain identity (3.1): T{x, y) ■T{y, z) = T{x, z), 

where the function T{., .) of two variables is defined as T{x, y) / + (x — y)K{x, y), 
or 

T{x,y)''^^ I + {x - y)F{xI - A)-^S-^{yI - B)-^G, ='(c\a^) x (c\c7b). 

Thus, we proved the following 

THEOREM 3.5. Let F, G be k x n and n x k matrices, A, B be n x n matrices with spectra 
a A SLnd as, and S be an invertible n x n matrix for which the equahty BS — SA = GF holds. 

Then the matrix function T(., .), which is defined by 

T{x, y) I+{x- y)F{xI - A)-^S-\yI - B)-'G, Vt = (c \ a a) x (c \ ae) , (3.28) 
satisfies the chain identity (3.1); 

T{x, y) ■ T{y, z) = T{x, z) 
and the diagonal unity identity (3.2): 

T{x,x) = I, 

and hence^^, is of form 

T{x,y) = R{x)R-\y), 
where R{x) =^ T{x, oo), {R{y))-^ = T{oo,y): 

R{x) = I - F{xl - A)-^S-^G, R-\y) = I + FS-\yI - B)-^G. (3.29) 

arc mutually inverse [i.e. R{x)R~^{x) = R~^{x)R{x) = I) rational matrix functions. 

PROOF. The chain identity (3.1) for the function T, defined by (3.28), was proved immediately 
before. The diagonal unity identity (3.2) evidently follows from the expression (3.28). The equal- 
ity T{x,y) = r(a;, oo) • T(oo,y)) is the special case of the chain identity (3.1) (written for the 
triple of points x,oo,y)). See Theorem 3.1. 

Letting y tend to oo, we obtain the expression (3.29) for the matrix function R{x) =^ T{x, oo). 
Letting x tend to oo, we obtain the expression (3.29) for the matrix function R~^ =^ T{oo,y). 

According to Theorem 3.1: the set (C \ (ta) H (C \ as) = C \ ((Ja U as) is not only nonempty but also very 
rich. So, we have many possibilities for choice of a distinguished point. ( See the proof of the Theorem 3.1). In 
particular, we can choose the point oo as a distinguished point. 
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That the functions i? and defined by (3.29)), are mutually inverse follows from the chain 
identity written for the triples x, oo, x and oo, x, oo. That the function R is rational is evident, g 

However, Theorem 3.5 says nothing about the nature of the rational function R. Imposing 
restrictions on the data A, B, F, G, we can say more about the matrix function function R. 

THEOREM 3.6. Let F be k x n matrix and G be n x k matrix with non-zero columns and 
non-zero rows respectively, i.e. 



G 
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where no column fi, f2, ■ ■ ■ , fn ^.nd no row gi, g2, • • • ,gn are zero , and let A, B be n x n be 
diagonal matrices with simple disjoint spectra, i.e. 

A = diag(Ai, ... ,A„), B = diag(/xi, ... ,)U„), 

where Ai, ... , A„; fj,i, ... , /i„ are pairwise different complex numbers. 
Assume that the n x n matrix S, 



S — ||'Spjg||l<pjg<n) 



9pfq 



(3.30) 



(which can be obtained from the data F, G, A, B as the unique solution of the Sylvester- 

Lyapunov equation BX — XA = GF ) is invertible. 

Then: 

1. The matrix function T{. , . ) of two variables, which is defined by 

T{x,y) = I+{x-y)F{xI-A)-^S-\yI-B)-'G, Vt= (c\(7a) x (c\(7b), (3.28) 
satisfies the chain identity (3.1); 

T{x,y)-T{y,z)=T{x,z), 
and the diagonal unity identity (3.2); 

T{x,x) = I. 

2. The matrix function T{x,y) is of the form 

T{x,y) = R{x)R-\y), 

where the matrix functions R, R"^ are defined by the formulas R{x) T{x,oc), R~^{y) =^ 
T{oo,y): 

R{x) = I - F{xl - A)-^S-^G, R-\y) = I + FS-\yI - B)-^G. (3.29) 
and are mutually inverse, {i.e. R{x)R~^{x) = R^^{x)R{x) = I). 

3. The matrix functions R and R~^ are rational matrix functions in general position. 
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4. The pole set V{R) of the function R coincides with the set {Ai, . . . A„}; the zero set J\f{R) 
of the function R coincides with the set {f^l, . . . i.e. 

A-p = A, Aj^ = B, (3.31) 

where A-p and A are the pole and zero matrices of the matrix function R. 

5. The semi-residual matrices Fp, Gp, Fj^, Gj\f of the matrix function R can he expressed 
in terms of the data F, G and of the matrix S (which in its turn is expressible from the 
data F, G, A, B): 

Fp = F; Gm = G- Fm = F-S-'^- Gp = -S'^ ■ G; (3.32) 



6. The right zero-pole coupling matrix S'^ and the left zero-pole coupling matrix for the 
matrix function R can be expressed in terms of the matrix S: 

S"- = S; = S-\ (3.33) 



PROOF. 

• Items 1 and 2 of Theorem 3.6 are already proved. (See Theorem 3.5). 

• Let's prove that the matrix functions R and R~^ are in general positions and investigate their 
singularities. The expression (3.29) for R may be written in the form 

l<j<n ^ 

with the matrix R\ . is of the form 

^A, = fj-Vj, 

where Vj is k-th row of the n x k matrix V '= —S~^G. From this expression it foUows that 
R is holomorphic outside the points {Ai, . . . , A„} and its inverse R~^ is holomorphic outside 
the points {fii, . . . ,//„}. Let's focus on the point Xj. There are two possibihties: or Vj = 0, or 
Vj 7^ 0. If Vj = then R\. = and hence, the function R is holomorphic at the point Xj. If 
Vj 7^ 0. than the matrix Rx. is non-zero, and has rank one. (We recall that, according to the 
assumptions of Theorem 3.6, fj ^ 0.) We show now that the equality Vj = is impossible. This 
equality may be written in the form EjS~^G = 0, where Ej = diag((5ij, 52j, ■ ■ ■ ,(>nj), { S is 
the Kronecker symbol). Multiplying the identity BS — SA = GF by the matrix EjS^^ from 
the left and by the matrix from the right and taking into account that the matrices ej and 
A are permutable (both are diagonal), we come to the equality {EjS~^)B — A{EjS~^) = 0. 
Because the spectra of A and B are disjoint, we obtain that EjS~^ = 0, and hence Ej = 0. The 
contradiction shows that the equality Vj = is impossible. Thus, each point Xj, j = 1, 2, ... ,n, 
is simple pole of the matrix function R, with residue matrices of rank one. 

• Analogously, we can show that the matrix function R~^ is holomorphic outside the points 
{Hi, . . . , Hn} and that each point fij, j = 1, 2, . . . , n, is a simple pole of the matrix function 

with residue matrices of rank one. Thus, items 3 and 4 of the claim of Theorem 3.6 are 
proved. 

• Item 5 of the claim follows from the representations (3.29). (Compare (3.29) with (1.37)- 
(1.38)). 



^^More precisely, one of the representatives of the equivalence class of the set of semi-residual matrices of the 
matrix function R. See Remark 1.2. 

'^^Which, in fact, serves as the definition of the matrix S. 
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• Now that we have estabhshed the relation (3.32) we may rewrite the equahty BS — SA = GF 
in the form Aj^S — SA-p = Gj\fF-p. Comparing the last equality with equality (2.52), we conclude 
that S'' = S. ^From (2.57) it follows now that 5' = 5"^ □ 

The "hybrid" version of this theorem (sec Remark 2.9 and formulas (2.66), (2.67)) can be for- 
mulated as well. This is the form which is convenient for applications in study of the Schlesinger 
system. 

THEOREM 3.7. Let F be k x n matrix and G be n x k matrix with non-zero columns and 
non-zero rows respectively, i.e. 



F=[hh---U 



G 
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where no column /i, /2, • • • , /n ^.nd no row gi, g2, . . . , Qn are zero , and let A, B be n x n be 
diagonal matrices with simple disjoint spectra, i.e 

^ = diag(Ai, ... ,A„), S = diag(;Ui, ... ,//„), 

where Ai, ... , A^; Mi, ... , fin are pairwise different complex numbers. 
Assume that the nxn matrix S, 



\^p,q\\'i-<p,q<ni 



_ 9pfq 



(3.34) 



[which can be obtained from the data F, G, A, B as the unique solution of the Sylvester- 

Lyapunov equation AX — XB = GF) is invertible. 

Then: 

1. Tiie matrix function T of two variables, which is defined by the formula 

T{x, y) = I-{x- y)FS-\xI - A)-^S{yI - B)-^S-^G, (3.35) 
satisfies the chain identity (3.1): 

T{x,y)-T{y,z) = T{x,z), 
and the diagonal unity identity (3.2); 

T(x, x) = I. 

2. The matrix function T{x, y) is of the form 

T{x,y) = R{x)R-'\y), 

where the matrix functions R, R~^ are defined by the formulas R{x) T{x,oc), R~^{y) =^ 
r(oo,2/): 

R{x) = I + FS-^{xI-A)-^G, R-\y) = I - F{yl - B)-^S-^G. (3.29) 
and are mutually inverse, {i.e. R{x)R~^{x) = R~^{x)R{x) = I). 
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3. The matrix functions R and R are rational matrix functions in general position. 

4. The pole set V{R) of the function R coincides with the set {Ai, . . . A„}; the zero set J\f{R) 
of the function R coincides with the set {fj,i, . . . i.e. 

A-p = A, Aj^ = B, (3.36) 

where A-p and Aj\f are the pole and zero matrices of the matrix function R. 

5. The semi-residual matrices Fp, Gp, Fj^f, Gj^j- of the matrix function R can he expressed 
in terms of the data F, G and of the matrix S (which in its turn is expressible from the 
data F, GAB): 

Fm = F; Gp = G; Fp = F ■ S-^; Gm = -S-^-G; (3.37) 

6. The right zero-pole coupling matrix and the left zero-pole coupling matrix for the 
matrix function R can be expressed in terms of the matrix S: 

= S; 5^ = S-\ (3.38) 

PROOF. Theorem 3.7 is nothing more then Theorem 3.6 in "other variables". Let A, B, F, G 
be the data of Theorem 3.7, and S be the matrix (3.34) generated by this data. Let's introduce 
the matrices 

F = FS-^, G = -S-^G, S = S-^. 

Equahty AS - SB = GF rewritten in terms of A,B, F, G, S becomes the form BS-SA = GF. 
^From the last equality it is easy to see that no column of the matrix F and no row of the 
matrix G equals zero: the equality EjG = or FEj = 0, where Ej = diag((5ij, 52j, ■ ■ ■ ,Snj), ( 
S is the Kronecker symbol), leads to the equality EjS = or SEj = 0, what contradicts to the 
invertibility of S. Now Theorem 3.6, applied to the matrix function 

T{x, y) = I+{x- y)F{xI - A)-\S)-\yI - By'G, 

gives the chain and diagonal unity identities for this T and the factorization T{x,y) = 
R{x)R~^{y), as well as the expressions for the semi-residual matrices Fp , Gp , Fj\f , Gj\f of the 
matrix function R: 

Fp = F, Gm = G, F_^ = F-S-\ Gp = -S'^G. 
Rewritten in terms of F, G, S, these relations becomes the form (3.37). □ 

4. THE SYSTEM REPRESENTATION AS A TOOL FOR THE SPECTRAL 
(WIENER-HOPF) FACTORIZATION OF MATRIX FUNCTIONS. 

In this section we show that the system representation may be used as an efficient tool for 
the so called spectral factorization (or the Wiener-Hopf factorization) of a matrix function. The 
problem of the spectral factorization can be formulated in the following way. 

GEOMETRIC CONFIGURATION. In the extended complex plane C a simple closed contour 
T is given. This contour separates C inter two regions, G+ and G- . These regions G+ and G- 
are connected open sets. We assume that the point {oo} does not belong to the contour T, thus 
one of the components, say G-, contains the point {oo}. 

^''See the footnote 
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DEFINITION 4.1. Given a k x k matrix function $ on the contour T, the factorization of the 
form 

4>(C) = 'J>+(C)-«'-(C) (Cer), (4.i) 

where and are kx k matrix functions, the matrix function $_|_ and its inverse ($+)~^ are 
holomorphic on G- U F, and the matrix function and its inverse ($_)~^ are holomorphic^^ 
on G+ U r, is said to be the spectral factorization (or the Wiener-Hopf factorization) of the matrix- 
function $ with respect to T. g 

We impose the normalizing condition 

$+(oo) = /. (4.2) 

on the factor (The function is holomorphic and invertible at the point oo, so the condition 
(4.2) makes sense). 

Under the normalizing condition (4.2), the spectral factorization (4.1) is unique. 

Even in the scalar case A; = 1 (i.e. $ is a complex valued function) the factorization problem 
(4.1) is not always solvable: there is a topological obstacle for the solvability. For a smooth 
nonvanishing complex valued function on F, the factorization problem (4.1) solvable if and 
only if there exists an univalued continuous branch of the function ln$(^) on F. If this condition 
is fulfilled, the solution of the factorization problem may be expressed in terms of the data $ by 
the formula 

*±(C) = exp{±^/i^*} (,eG±). (4.3) 

r 

The proof of the fact, that the formula (4.3) gives the solution of the factorization problem 
(4.1) in the scalar case, is based essentially on the Sokhotskii-Plemelj formulas on the boundary 
behavior of the Cauchy integral. Actually, in the scalar case we solve the additive problem 

*+(C) + *-(C) = ln$(C) (CeF), 

and then we exponentiate. In the matricial case we still can solve the appropriate additive 
problem using the Cauchy integral, but exponentiating does not lead to the desirable result: 
In view of noncommutativity of the matricial multiplication, exp{^ + B} ^ exp{A} ■ exp{B} 
for matrices A and B in general. In the matricial case, the situation with the factorization 
problem (4.1) is much more complicated than in the scalar case. There are not only topological 
obstacles to the solvability of this problem. The factorization problem (4.1) is equivalent to 
some system of singular integral equations on F, and to analyze this system is approximately so 
hard as to investigate the original factorization problem (4.1). The factorization problem (4.1) 
appeared firstly in the context of Hilbert's twenty-first problem: to construct a Fuchsian linear 
differential system with the prescribed monodromy group. See [P12], [Bo] and [Gah] for details 
and historical references. The factorization of the type (4.1) is used also in the solving systems of 
singular integrals equations with Cauchy kernel on the contour as well as for solving of systems 
of integral equations which kernel depends on the differences of the arguments on the half axis. 
Sec [Vck], [GoKr] and [ClGo] on this subject. It should be mentioned that the pioneer papers 
[Birl] and [Bir2] by G. Birkgoff had a profound impact on the further investigations on matrix 
factorization. 

We confine ourself to the case the function $ is a rational matrix function (or, more precisely, 
the restriction on F of a rational matrix function) such that the functions <5 and are 

'^^In particular, the functions $+ and $_ are holomorphic on the common boundary T of the domains G+ and 
G-, so the relation (4.1) makes sense. 
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holomorphic on the contour T. In this case the factorization (4.1) is global, i.e. the matrix 
functions and are rational, and the equahty 

= $+(z) • (Vz G C) (4.4) 

holds. Indeed, in this case the function ^.j.^'^^ is holomorphic within G_ except finite many 
poles located on the set P($) n G-. In view view of (4.1), this function continues analytically 
into the function which is holomorphic on G+ U T. Thus, the matrix-function <!>_ has no 
other singularities in C than finite many poles and hence is rational. For the same reasoning, 
the matrix function is rational. 

Thus, in the case that the initial matrix function $ is rational, the problem of the spectral 
factorization may be reformulated in the following manner: 

DEFINITION 4.1'. Given a k x k rational matrix function <I>, det<& ^ 0, its factorization of 
the form (4.4), where are rational matrix functions with zero and pole location 

P($+) c G+, Afi^+) C G+, P($_) c G_, A/'($_) c G_, (4.5) 

is said to he the spectral factorization (or the Wiener-Hopf factorization) of the matrix-function $ 
with respect to F. □ 

We consider even the more special case: the function $ is a rational function in general position. 
In this case the calculation of the factors and can be performed completely by hand, in 
terms of poles and "zeros" of the matrix function $ and its semiresidual vectors. 

So, let $ be a rational matrix function in general position, normalized by the condition 

$(oo) = /. (4.6) 

Let 7^($) and 7V(<I>) be its pole and zero sets, j4p($) and Aj^{^) be its pole and zero matrices, 
F'p{^), Gp(<I>), F_A/($), Gj\f{^) be the appropriate semiresidual matrices. According to the The- 
orem 2.6, the zero-pole coupling matrices S^'{^) and S^{^) are invertible, and the matrix-function 
$ admits the representations of the form (2.60) and (2.62): 

$(z) = I - Fp{^){zl - At,^^))-^S^{^)-'Gm{^) (4.7) 

and 

$(z) = / + Fv($)5'($)~\z7- Ap($))-^Gp($). (4.8) 

The Sylvester-Lyapunov equations for the matrices S'^{^) and S\^) (which actual-y are the 
definitions of these matrices) are of the form: 

A^($)5''($) - 5'"($)Ap($) = G^($)Fp($), (4.9) 

and 

Ap($)5'($) - S\^)A^f{^) = G'p{^)F^r{^). (4.10) 
Moreover, the matrices S"'(^>) and S''(<I>) satisfy the equality 

S\^)-S''{^) = I. (4.11) 

(This is (2.57) for the matrix function $). The zero-pole coupling relation hold: 

a). Gaa($) = -S'{^)Gv{^); b). Gp($) = -S\^)Gm{^); 

(4.12) 

c). Fp($) = FAr(^>)5"($); d). Fj^{<^) = Fr{^)S\^). 
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(This is (2.51) for the matrix- function 
According to the assumptions, 

a). -P($)nr = 0; b). Ar($)nr = 0. (4.13) 

Assume that the factorization (4.4) holds, and that the normalizing condition (4.2) is satisfied. 

Since the functions are holomorphic in and the functions $, have only 

simple poles in G+, from the relations 

a). $+ = $-<^>l\ b). $-i = $_-$-i (4.14) 

it follows that the functions <I>_|-, <1>^^ have only simple poles in G+. In G_ and on T the functions 
and don't have singularities at all. Thus, 

P($+) =P($)nG+, Mi^+)=M{^)nG+. (4.15) 

Let A e V{^+). From (4. 14. a) it follows that the residues R\{^+) and R\{^) of the matrix 
functions and $ at the point A are related by the equality 

i?A($+) = i?A(^)-(^-(A))-^ (VAgP(<&+)). (4.16) 

Since $ is a matrix function in general position, the rank of the residue matrix Rx{^) is equal 
to one. According to the assumptions, the matrix ^*-(A) is invertible for A G G+. ^From (4.16) 
it follows now that the rank of the residue matrix Rx{^+) is equal to one as well. Let now 
/Li G A^($+). /.From (4.14.b) it follows that the residues i?^($+) and Ri^i^) of the matrix 
functions and at the point /i are related by the equality 

i?^($+) = $_(/x) • i?^($) ( V/x G Ar($+) ). (4.17) 

The rank of the matrix R^i^) is equal to one ($ is a matrix function in general position); the 

matrix ^*-(/i) is invertible (according to the assumptions, the matrix ^*-(/i) is invertible for 
jj, G G+). From (4.17) it follows now, that the the rank of the residue matrix R^^^^) is equal 
to one as well. ^From (4.15) it follows that 

P($+)nAr($+) = 0; {oo}^r{^+), {oo}^ J\f{^+). 

Thus, is a rational matrix function in general position. 

In the same way we obtain that is a rational matrix function in general position, and 

P($_) =P($) nG_, Af{^-) =Af{^)nG-. (4.18) 

Moreover, the residues Rx{^-) and i?^($_) of the matrix functions and $1^ at the poles 
A G r{^-) and G A/'($-) ( = V{^Z^) ) are related to the residues of the matrix functions $ 
and at the same points by the equalities 

Rx{^-) = ($+(A))-^ • Rxm ( VA G P($-) ) (4.19) 

and 

i?^($_) = i?^($) • $+(/x) ( V/x G M{^+). (4.20) 

^From (4.16) and (4.17) it follows that the left semiresidual vectors of the matrix functions 
and $ at the poles A G coincide: 

/a(^>+) = /a(^>) (VAGn^+)). (4.21) 
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and the right semiresidual vectors of the matrix functions and $ ^ at the poles fi G 
coincide: 

g^{^-^) = g^{^-^) (V^GAA(cI>+)). (4.22) 

In the same way we can obtain that the right semiresidual vectors of the matrix functions $_ 
and $ coincide: 

gx{^-) = 9x{^) ( VA G ). (4.23) 

and the left semiresidual vectors of the matrix functions $1^ and coincide: 

M^Z') = U^-') (V//GA/-($_)). (4.24) 

The eqialities (4.21) - (4.24) are crucial for solving of the considered factorization problem. 

According to Lemma 1.2, the equality 

#P($+) = #Ar($+) 

holds for the rational matrix function in general position $_|_. Taking into account the equality 
(4.15), we obtain the following equality 



#(P($) n G+) = #{M{^) n G+) n+. 



Of course, the equality 



#(P($)nG_) = #(Ar($)nG_) 



def 



ri-. 



(4.25) 



(4.26) 



holds as well. 

To simplify notations, we denote 

r+ = P($) n G+ ; V- = P($) n G_ ; Af+ = AA($) n G^ 



AA_ = A/'($) n G_ 



To the decompositions 



p($) = v+uV-, v+nV- 



(4.27) 



of the pole and zero sets 'P{^), AA($) of the matrix function $ there correspond natural 
block-decompositions of the matrices which appear in the system representations of the ma- 
trix functions $ and $~^: pole and zero matrices A-p[^), Aj\f{^), the semiresidual matrices 
F'p{^),Fj^{^),G-p{^),Gj^{^) as well as the zero-pole coupling matrices 5"'($) and 5'($). 

Namely, the decompositions of the pole and zero matrices ^•p($) = diag(A;);yj£-p($) and 
Ajv'($) = diag(/x;)^jgjv'(*) (of the dimension n x n) are of the form: 



A-pi^)! 
A-p{^)2 



(4.28) 



A^f{^)l 
Aa/-($)2 



(4.29) 



where 



A-p{^)i = diag{\i)xiev+; A-p{^)2 = diag(Ai)A,eP_; 
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Am{^)i = diag(/x;)^,6jv-+; Aj^{^)2 = diag(w)^,ejv--; 

A-p(^)i and Aj^(^)i are diagonal n+ x matrices, A'p{^)2 and ^a/'($)2 are diagonal n_ x n_ 
matrices; and n_ arc defined in (4.25) and (4.26): n+ = #7^+ = t^A/+; n_ = = ^M-). 

The block-decompositions of the semiresidual matrices Fp($), Fv'($), G-p{^), Gjv'(^): 

Fp($) = row(/AjA,GP(<i')> ^JV(*) = row(/^J^je^($) are A; x n matrices, 
Gp($) = col(5rAjA,GP($), Gaa($) = col(5(^J^,e^($) are nxk matrices. 



related to the matrix- function are of the form: 



p). Fri^) = [Fp($)i F-pi^)2] 



n). Fa/-($) = [Faa($)i Fj^{^)2] 



and 



p). Gp(#) 



Gp(<i>)i 



n . 



Ga/-($)i 
GAr(^')2 



(4.30) 



(4.31) 



where 



F-p($)i = row(/AjA,eP+, ^V(^)i = row(/^,)^;e^^ are k x n+ matrices; 

F'p{<^)2 = row{fxi)xieV-, -^Ar(^*)2 = row(/^J^;e^_ are k x n_ matrices; 

G-p($)i = col{g\^)xier+, G^f{^)l = col(5(^,)^;eArf are n+ x k matrices; 

0^(^)2 = col{gxi)xieV-, Gjv-(^)2 = col(c/^,)^,ejv-- are n- x k matrices. 



iProm (4.15) and (4.18) it follows that 

p). ^p($+) = ^p($)i 

and 

p). Ap{^_) = A-p{^)2 
In view of (4.21) and (4.22), 

p). Fp($+) = Fp($)i, 

In view of (4.23) and (4.24), 

p). Fm{<P_) = Fm{^)2, 



n). A^{^+) 
n). Aj^{^_) : 



Aj^{^)2- 



n). G^($+) = G^($)i. 



n). G7,($_) = Gp($)2. 



(4.32) 
(4.33) 

(4.34) 
(4.35) 



Thus, we have expressed the pole matrix A'p{^^) and semiresidual matrices F-p($+), Gj\/{^^) 
for the left factor in terms of blocks of the appropriate block-decompositions of the pole 
matrix ylp($) and and semiresidual matrices F-p{^), G_xf{^) for the factorized matrix function 

We also have expressed the zero matrix Aj\f{^-) and semiresidual matrices (<!>_), G-p($_) 
for the right factor in terms of blocks of the matrices Aj\/{^), Fj^{^), G'p{^). In principle, 
these data are sufficient to recover the factors and ( from the appropriate blocks of the 
pole and semiresidual matrices for the factored matrix function $ ). To carry out the recovering, 
we have to solve the Sylvester-Lyapunov equations (2.52) and (2.53) for the matrix functions 
and respectively to find from these equations the zero-pole coupling matrices S*" (<!'+) 
and S\^-): 



AAr(<l>+)5"($4 
Ap{^_)S\^. 



5'-($+)^,.($+) 

-5'($_)y4p($_) 



GAr($+)Fp($H 
Gp($_)FAr($- 



(4.36) 
(4.37) 
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According to the assertion 1 of Theorem 2.6 (actually, according to Theorem 2.5: see (2.31)), 
their solutions S^{^+) and S'^($_) are invertible matrices. Then we construct the factors $_|_ and 
$_ according to the formulas (2.60) and (2.62): 



$+(^) = / - F-pi^+){zI - A-p{^+))-\S'-i^+))-'Gx{^+) 



(4.38) 



and 



$_(z) =1 + Fm{'^_){S^{'^_))-\zI - Ar{'^-)y^Gp{<^-). (4.39) 

We express now the matrices S^{^+) and 5'($_) in terms of blocks of the matrices S^{^) and 
S\^). The block-decomposition of the zero-pole coupling matrices S^{^) and S\^) are of the 
form: 

5^($)21 



5''($) 



(4.40) 



(4.41) 



where 11- block-entries S^{^)u, S''($)ii are n+ x n+ matrices, and 22- block-entries S"'($)22, 
'S''('I*)22 arc n_ x n_ matrices (n+, n_ arc defined in (4.25) and (4.26)). The block- 
decompositions (4.40) and (4.41) are consistent with the block-decompositions (4.28)-(4.31) 
of the pole and semiresidual matrices. The Sylvester-Lyapunov equation (4.9) for the matrix 
S'''($), written in the block-matricial form 



Ax{^)i 
A^{<^)2 



5'-(1>)ii 

S'{^)21 



5'-($)22 



5'-(a>)n 
5'-($)2i 



5"(1>)l2 
^'■(1>)22 



ylp($)i 






Av{^)2 



Gm{^)2 



[Fp($)i Fp($)2], 



may be considered as a system of matricial equations for the block-entries of the matrix S'^'{^). 
This system is decomposed into four equations for block-entries of this matrix. 

In particular, the equation for the entry S^{^)ii is of the form 



Am{<^>) iS'i'^)!! - S'i<i>)iiAvi<i>) 1 = GAr(1>)iFp($)i. 



(4.42) 



The Sylvester-Lyapunov equation (4.10) for the matrix S''($), written in the block-matricial 
form, may be considered as a system of matricial equations for the block-entries of the matrix 
S''($). In particular, the equation for the entry S''($)22 is of the form 



Ar{<^)2S\<^)22 - 5'(^>)22^^(^)2 = Gp($)2F^($)2. 



(4.43) 



We show that 
and 



5'-($+) = 5'-($)n (4.44) 

S\^-) = S\^)22. (4.45) 

The easiest way to do this is to use the explicit formulas (2.56) for the matrices S'', S"' in terms 
of the poles, zeros and semiresidual matrices of the matrix function. For the function $ the 
formula (2.56.r) takes the form 



5/.p($)-/A,($) 



(4.46) 
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and 



In particular, 

S'{^)n = K,,m,,eMi^)nG^,x.evmnG^, ^;,($) = ^^^^^^ ' ^'■^^^\ (4.47) 
For the function the formula (2.56.r) takes the form 

= ||4,,($+)IU,.AA(*+),A,.P(*+), = ^^-^^f'Jl"^^'-^ (4.48) 

jjjp Aq 

Comparing two last formulas and taking into account (4.15), (4.21) and (4.22), we conclude that 
(4.44) holds. In the same way, comparing the formulas 

Ap fJ,q 

5'($_) = ii4,,($-)iiA,eP(*-),^.6AA(*_), = ^"''^^r^/r^'^"^ (^^-^o) 

and taking into account (4.18), (4.23) and (4.24), we conclude that (4.45) holds. 

The explicit expressions (4.48), (4.50) for the zero-pole coupling matrices are consequence 
of the Sylvester-Lyapunov equalities (4.36), (4.37). It is also possible to derive the equations 
(4.44) and (4.45) directly from the equalities (4.36), (4.37), bypassing the explicit expressions 
(4.48), (4.50). (4.46), (4.49). The latter way is better, because it is applicable not only to matrix 
functions in general position. 

Taking into account the equalities (4.32), and (4.34), we came from (4.36) to the equation 

AAr($)i5'-($+) - 5'-($+)^p(#)i = GAr($)iFp($)i. (4.51) 

Comparing (4.42) and (4.51) and taking into account the uniqueness of the solution of the 
Sylvester-Lyapunov equation, we conclude that the equality (4.44) holds. 

In the same way, we can establish the equality (4.45). Taking into account the equalities 
(4.33), (4.35), we came from (4.37) to the equation 

^p($)2-S'($_) - S\^-)Am{^)2 = Gv{^)2Fm{^)2. (4.52) 

Comparing (4.43) and (4.52) and taking into account the uniqueness of the solution of the 
Sylvester-Lyapunov equation, we conclude that the equality (4.45) holds. 

According to Theorem 2.6 (applied to the rational matrix functions in general position 
and ^>-), the matrices 5^($+) and S''(^>_) are invertible. In view of (4.44) and (4.45), the block- 
entries and (5'(*))22 (of the matrices S""($) and S\^) respectively) are invertible. 

Now we may rewrite the formulas (4.38), (4.39) for the factors <I>-|_ and <!>_ in terms of block- 
entries of the pole-, zero-, semiresidual- and zero-pole coupling matrices for the factorized matrix 
function Substituting the expressions (4.32.p), (4.34) and (4.44) into (4.38), we obtain 

$+(z) = / - Fp($)i(zl - Ap($)i)-i(5'-($) ii)-^G^($)i. (4.53) 

Analogously, substituting the expressions (4.32. n), (4.35) and (4.45) into (4.39), we obtain 

^^{z) = I + Fu{^\{S\^)^^)-\zI - Ar{^\)-^Grm^. (4.54) 
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The formula (4.53) expresses the factor in terms of the values G_\f{^) and S'^{^), 

which appear in the representation (4.7), whereas the formula (4.54) expresses the factor $_ in 
terms of the values Fj\f{^), Gp($) and >S'^($), which appear in the other representation (4.8). 
This disagreement is inconvenient for some calculations. Therefore we also give the formula 
which expresses the factor $_ in terms of the the same values F^($), G^{^) and 5^($), which 
appear in the expression (4.53) for the factor To do it, we have to use the zero-pole coupling 
relations (4. 12. a), (4.12.c), Together with (4.35.p) and (4.34.n), these relations mean: 

Fv($_) = Gp($_) = -((>S'^($))-iGAr($))2. (4.55) 
$): The equality (4.11) together with (4.45) means: 

5'($_) = (S'-($)-i),2. (4.56) 
Substituting the expressions (4.55) and (4.56) into (4.39), we obtain 

$_(z)=7-(Fp($)5'^($)-i)^((5'^($)-i) J"'(z/-^p($)2)"'(5^($)-^Ga/-($))^. (4.57) 

The last formula already expresses the factor in terms of the values Fp($), Gj\f{^) and 
S"'($), i.e. in terms of the same values which appear in the representation (4.53) of the factor 

However, for the further considerations it will be useful to transform this formula, substituting 
into it the expression for the inverse matrix 5"'($)~^ in terms of the block-entries of the matrix 
S^{^) itself. First of all we recall a formula for the inversion of a 2 x 2 block matrix with square 
diagonal block-entries. Let 



M 



(4.58) 



mil rni2 

be a square n x n block-matrix matrix, which block-entries mn and m22 be square ni x ni and 
n2 X n2 matrices respectively (n = ni -|- 722). We assume that the matrix mu is invertible. To 
inverse the matrix M, we first of all factorize it: 



mil 


mi2 




_ m2i 


m22 





The matrices 



m2imii 
I ' 



m2i m^/ / 




I 



and 



mil 








,-1 , 



m^i mi2 







m22 - m2i m^]^ mi2 

are invertible, and 



(4.59) 



mif^mi2 



/ 

m2imri^ 




/ 





I 



I m,ylmi2 
/ 



-m2\m^^ 

Therefore the matrix M is invertible if and only if the matrix 







mil 



-mi/mi2 







-1 , 



m22 — m2i m-^{ mi2 
-1 , 



IS 



invertible. The latter is invertible if and only if the matrix (m22 — m2i m^^^ mi2) is invertible. 
If the matrix (m22 — m2i rn^^ mi2) is invertible then 







-1 


mil 


mi2 




_ m2i 


m22 





/ —m-^lmi2 



m 



11 







(m22 — m2i m^/ mi2 



-m2i m-^l I 



The invertibility of the matrix mii was assumed from the very beginning 
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Finally, for the matrix M with the decomposition (4.58), 



In particular. 



m^"/ 




+ 



I 



m22 — m2i m-^imi2 



—17121 m^x I 



(M ^)22 = ("1-22 - m.2i m^/ mi2) S 



(4.60) 



(4.61) 



\i F \s k X n matrix, with the block-decomposition F = [Fi F2], where Fi, F2 are k x m and 
k X n2 matrices respectively, then 2-entry of the matrix F M^^ is of the form 



(FM-i)2 = [Fi F2] 



-m^^mi2 
I 



■ ("7,22 — "721 m-i/ "1-12) ^ ■ 



(4.62) 



If G is n X A; matrix, with the block-decomposition G = 



Gi 
G2 



where Gi, G2 are ni x k and 



772 X A; matrices respectively, then 2-entry of the matrix M ^ G is of the form 
[M~^G) 2 = (77722 — "721 ruil 77712) ^ • [-fn2im]^l l] ■ 



Gi 
G2 



(4.63) 



Let us take the matrix S'^($) with the decomposition (4.40) as the matrix M {n\ = 77+, 772 = 77_) 
as the matrix Af. Both the matrix S"'($) and its block-entry the matrix S'^($)ii are invertible. 
Therefore the entry (S'''($)~^)22 of the inverse matrix {S^{<^))~^ = S\^) is invertible as well, 
and, according to (4.61), 



{{s'^m-')22)-' = s'm22 - s'm2i 5'-($)r/ s'mi2. 

In particular, the matrix Sf^(*)22 - 'S'^(*)2i S^'i^)^! 5f^(*)i2 is invertible. 

In view of (4.11), the equality (4.45) can be presented in the form 

5^($_)22 = (5^(^)22 - S'\^)21 5^($)2-/ S^muV' . 



(4.64) 



(4.65) 



Taking the matrix Fp(<I>) with the decomposition (4.30.p) as the matrix F and the matrix 
Gj\f{^) with the decomposition (4.31. n) as the matrix G. we reduce the formula (4.57) to the 
form 



$_(z) = 7-[FAr($)i F^m 



I 



■ (z7-^p($)2) 



■{s'\^)22 - S^i^)2i 5''(^)n '5''(^)i2) • [-S'\^)2i S' ^n^ l\ 



G^mi 
Gjv($)2 



(4.66) 



Thus, we proved the following 



THEOREM 4.1. Let $ be a rational k x k matrix function in general position, satisfying the 
normalizing condition (4.6); <I>(oo) = /, with pole and zero sets 'P(<1>) and 7V($), pole and zero 
matrices A-p{^) and Aj\/{^), semiresidual matrices Fp{^), Fjv'($), G-p($), Gj\f{^), and zero- 
pole coupling matrices S"'($) and S'^($), which are decomposed into blocks as described above. 
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Assume that the function $ admits the spectral factorization with respect to T, i.e. 

$(z) = $+(^) (zeC), (4.4) 

where are rational matrix functions, with the pole- and zero-location: 

V{^+) c G+, c G+, c G_, c G_. (4.5) 

Let the normalizing condition (4.2); $_|_(oo) = I hold. 
Then: 

1. For the matrix function $, the number of poles and the number of "zeros", located in 
G+, are equal: #(P(<I>) n G+) = #(A/'(<I>) n G+); the number of poles and the number of 
"zeros", located in G-, are also equah #(P($) n G_) = fl G_). 

2. The block-entries S^{^)n and /S''($)22 {in the above described block decompositions 
(4.40), (4.41) of the matrices S^{^) and S\^) respectively), as well as the matrix 
S'{^)22 - S'{^)2iS'{^)i^S'{^)u axe square invertible matrices. 

3. The factors and are rational matrix functions in general position, which are repre- 
sentable in the form (4.53) and (4.54) respectively. 

4. The representation (4.53) is the right system representation of the factor the entries of 
this representation is expressible in terms of the entries of the right system representation 
for (4.15), (4.32), (4.34) and (4.44) hold. The representation (4.54) is the left system 
representation of the factor the entries of this representation is expressible in terms of 
the entries of the right system representation for (4.18), (4.33), (4.35) and (4.45) hold. 
The representation (4.54) may be rewritten in terms of the entries of the right system 
representation for $_|_; (4.66) holds. 

The converse statement is true as well. 

THEOREM 4.2. Let $ be a rational k x k matrix function in general position, satisfying 
the normalizing condition (4.6): $(oo) = /, with pole and zero sets 'P(<I>) and A^($), pole and 
zero matrices A'p{^) and A^{^), semiresidual matrices F'p{^), Fj\f{^), Gp(<I>), Gj\f{^), and 
zero-pole coupling matrices S^{^) and >S'("3>), which are decomposed into blocks as described 
above. 

Assume that the following conditions are satisfied: 

1. JVone of the poles of ^ as well as none of its "zeros" belongs to T: 

p($)nr = 0, p($-^)nr = 0. 

The number of poles and the number of zeros of the matrix-function located in G+, are 
equal: 

#(P($)nG+) = #(A^($)nG+), 

or, what is the same, the number of poles and the number of zeros of the matrix-function 
$, located in G_, are equal: 

#(P($)nG_) = #(Ar($)nG_). 
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2. The block-entry of the matrix {in the above described block decompositions 

(4.40)) is invertible matrix. 

Then: 

1. The matrix functions and which arc defined by (4.53) and"^^ (4.54), are rational 
matrix functions in general position. For these functions and the conditions (4.5) 
are satisRed. 

2. The matrix function $ admits the spectral factorization (4.4) with respect to F, with these 
matrix-functions and as the factors. 

PROOF. To investigate the properties of the matrix functions $_)_, we will use Theorem 3.6. 
To this theorem be applicable to the functions we have to check whether its assumptions 

follow from the assumptions of Theorem 4.2. First of all, we have to be sure that none of the 
columns of the matrices F'p{^)i, -Fa/(<I>)2 is a zero-column and none of the rows of the matrices 
Gj\f{^)i, G-p($)2 is a zero-row. This property holds indeed because it holds for the including 
matrices F-p($), Fv($), G-p($), Gj\f{^) (as for semiresidual matrices for the rational matrix 
function $ in general position). Then we have to check that the matrices S"'($)ii, which appears 
at the "core" of the representation (4.53), satisfy some Sylvestcr-Lyapunov equality of the form 
BS''{^)n - 5"'($)ii A = Gj^{<^)i ■ Fp($)i, where A = yl-p($)i and 5 is a diagonal matrix. In 
the considered case, the identity (4.42) plays the role of such Sylvester-Lyapunov equality, with 
B = Aj\f(^)i. As we are already established, the equality (4.42) is the consequence of the equality 
(4.9). The latter holds as the Sylvester-Lyapunov identity for the matrix function Thus, the 
representation (4.53) of the function $_|_ is its system representation. According to Theorem 
3.6, the matrix-functions is a rational matrix-function in general position, which pole- and 
zero-sets are totalities of the diagonal entries of the matrices A-p{^)i and Ax{^)i respectively: 
V{^+) = jdiag A7r,($)i|, and Ar($+) = |diag^A/-($)i|. Hence (see the definition (4.28), (4.29) 
of the matrices A{V)i, A{M)i ), the conditions P($+) = P(#) n G+ Af{^+) = M{^) fl G+ hold. 

The matrix S"'(<I>) is invertible as the Sylvester-Lyapunov matrix related to the matrix- 
function $ (Theorem 2.5. See (4.11)). The corner block-entry 5'''($)ii is invertible by the 
assumptions of Theorem. Hence (see the identity (4.60) with S'^($) as M), the matrix 
S^{^)22 - 5^($)2i>S'^($)2/S''($)i2 is invertible as well, and (4.64) holds. In view of (4.11), 
(4.65) holds. In particular, the block-entry S''(<5)22 is invertible as well. Thus, the matrix- 
function $_ is well defined by the formula (4.54). The Sylvester-Lyapunov equality (4.43) for 
the matrix S'^($)22 follows from the equality (4.10). (The latter holds as the 'left" Sylvester- 
Lyapunov equality for $). According to Theorem 3.6 (to be more precise, according to the 
"left" version of this Theorem), is a rational matrix function in general position, which 
pole- and zero-sets are totalities of the diagonal entries of the matrices A'p{^)2 and Aj\f{^)2 
respectively: 7-'($_) = |diag Ap(<l>)2|, and A/'($-) = |diagyl;v'(^)2|- Hence, the conditions 
p($_) = 77($) n G- A/'($_) = J\f{<^) n G- hold. All the more, the conditions (4.5) hold. 

Thus, if the equality (4.4) holds for $_|_, defined by (4.53), (4.65), it gives the spectral 
factorization of 

One remains only to verify the equality (4.4). The representation (4.54) is convenient to 

investigate the properties of However, to verify the equality (4.4), it is more convenient to 
use the representation (4.66). To derive (4.66) from (4.54), we have to use the zero-pole coupling 
relations (4.12), the formula (4.60) for inversion of 2 x 2 block- matrix (applied to the matrix 
S'^($) with the block-decomposition (4.40)) and, in particular, the equality (4.65). 

^°If the block-entry S"^($)ii is invertible, the block-entry S'{^)22 is invertible as well. 
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To abbreviate the notation, we omit some notation entries, like indices etc. So, for ex- 
ample, we write S instead S""(^», 5ii instead S'^{^)ii, F = [Fi F2] instead (4.30.p), G = 



G2 



instaed (4.30.n),^ 



Ai 
A2 



instead (4.28), B 



Bi 
B2 



instead (4.29). 



Using the 2 x 2-block-matrix inversion rule (4.60), we present the representation (4.7) for $ 
in the form 



-52i5n I ]g, (4.67) 



where 

A = S22 — S21S11S12. 
The representation (4.53) for the function we present in the form 



■( 






^11 ^12 











I 





$+(z) = / - F{zl - A)-^ 







G. 



(4.68) 



The representation (4.54) for rewritten in the form (4.66), is: 



$_(z) = I-F 



hi -^12 
I 



1 A-l 



{zl - A2)-'A 



S21S11 I 



G. 



(4.69) 



Multiplying the expressions in the right hand sides of (4.68) and (4.69) term by term, we obtain: 
^+{z)^-{z) = I- 



-F{zl - A)-^ 

-S11S12 
I 

+F{zI - A)~^ 







{zI - A2)-^A 



1 A-l 



1 




G-F 


— S-j^^Si2 









I 



G- (II) 
-52i5iY i]g+ (III) 

{Zl - ^2)-'A-l [ -S2lS^i^ I ] . (IV) 

(4.70) 

Substituting the expression for GF from the Sylvester-Lyapunov identity GF = BS — SA, or 

GF = S{zl -A)- {zl - B)S, 
into the expression (IV), we obtain: 

(IV) = 

°1 - ^ -^^''^''\zI-A2)-^A-^\ -S2iS^,' I 



F{zl - A)-^ 
-F{zl - A)-^ 








S{zl - A) 
■ {zl - B)S ■ 



•Jii '-'12 
I 

—S11S12 
I 



1 A-l 



{zI-A2)-'A 



S21SII I 



(4.71) 
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It is clear that 







{zl - B)-' 



* 


, s 






' 







I 




* 



Thus, 



{zI-B)-^S 



Sn 


the second summand in the expression (4.71) for (IV) vanishes, and 



'S'II^'S'12 




" " 


/ 








(IV) = F{zl-A)-^ 







S{zl-A) 



—S11S12 
I 



1 A-l 



{zI-A2)-^A 



S21SII I 



G. 



(4.72) 



As 





•5 = 7- 


—S^iSi2 







I 



(IV) = F 



-S11S12 
I 



-F{zl - A)-^ 
It is clear that 



—S-^^Si2 
I 



{zl - A2)-^A-^ 

(zI-A) 



—82181^ I 



G- 



-8-1^1 S12 
I 



izI-A2)-'A-^ -52i5iY / 



G. 



{zl - A) 



-Si^Si2 

I 



izi-A2r' 



I 



hence 



Thus, 



—Si^Si2 
/ 



{zl - A) 



-Si^Si2 

I 



{zi-A2r' = 



-S-j^lSi2 

I 



(IV) = F 
-F{zl - A)-^ 



811 S12 
I 



{zI-A2)-^A 



1 A-l 



—8218-^1 I 



G- 



—8]^^8i2 
I 

iProm (4.70) and (4.73) it follows that 



(4.73) 



-5'2i5'ii^ I 



G. 



$+(z) $_(z) = 7 - F{zl - A)-^ 



( 


- a-l 
•^11 


" 


+ 


—8^1 812 


A 












I 





—8218-^^1 I 



G. 



Comparing the last expression to (4.67), we conclude that the factorization equality (4.4) 
holds. Theorem 4.2 is proved. □ 

REMARK 4.1. Let $ be a rational matrix function in general position, satisfying the normal- 
izing condition (4.6): $(cx)) = 7. Assume that none of the poles of $ and none of the "zeros" of 
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$ belongs to the contour F, and the conditions (4.25) : #(P($) n G+) = #(AA($) n G+) is sat- 
isfied; (or, what is the same, the condition (4.26): #(P(^>) nG_) = #(A/'($) nG_) is satisfied). 
This means that the matrices 5"'(^>)ii and S''($)22 are square one. According to Theorems 4.1 
and 4.2, the matrix function $ admits the spectral factorization with respect to F if and only if the 
matrix S'^{^)ii is invertible, or, what is the same, the matrix 8^^)22 is invertible. What can wc say 
if the condition of the invertibility of the matrix S^{^)ii is violated? According to the matrix 
factorization theory, the matrix function $ admits the factorization of the form 

$(z) = $+(z) D{z) ^-{z), (4.74) 

where the matrix functions and are holomorphic on the set G_ UT, the matrix functions 
and $1^ are holomorphic on the set G+ U T and D{z) is the matrix function of the form 

D{z) = diag ((z - zor\ {z - zor, . . . , (^ - ^0)"*), (4.75) 

Ki, K2, . . . , Kk arc integer numbers. The numbers ki, K2, ■ ■ ■ , Kk are said to be the partial indices 
of the matrix function $ with respect to the contour T. 

The point zq is an arbitrary chosen fixed point from G+. Of course, the factors 
and D from the factorization (4.74) depend of the choice of the point zq. However, the set 
{ki, K2, ■ ■ ■ , Kk} (and under the normalizing condition 

1^1 < K2 < ... Kk (4.76) 

the matrix D) is determined uniquely by the function In particular, the partial indices do 
not depend on the choice of the distinguished point zq. In contrast to this, even for the given 
distinguished point zq the matrices $_|_ and are determined non-uniquely. However, this 
non-uniqueneess can be easily described. (See Theorem 7.1 from [GoKr]). 

The partial indices play a fundamental role; in the homogeneous Hilbert problem they were 
first introduced by N.I.Muskhelishvili and N.P.Vckua [MuVe]. (The factorization of another form: 
$(z) = ^j^{z) ^^{z) D{z) was considered by G.Birkgoff much earlier, in 1913 (see [Birl]). The 
matrices and D have the same properties that in above described factorization (4.74) 

(in particular, the matrix D is of the form (4.75)), but the matrices themselves are different. 
In particular, Birkgoff's partial indices may be different. The relation between these kinds of 
factorizations is studed in [FM]. 

The natural question arise how to compute the factors <!>_, D in terms of the pole, zero and 
semiresidual matrices for the given rational matrix function $ in general position? Of course, there 
are general methods for performing such a factorization (see, for example, [ClGo]). However, 
we believe that for a rational matrix function in general position the factorization may be done 
much more explicitly. q 

REMARK 4.2. Let as assume now that the rational matrix function $ in general position 
depends holomorphically on some parameter a € C and satisfy the normalizing condition (4.6): 
$(00, a) = /. The factors D in the factorization (4.74): 

$(z,a) = ^+{z,a)D{z,a)^-{z,a), (4.74 «) 

D{z, {a)) = diag ((z - zor^^'l {z - zor^"'\ • • • , (z - zo)"'^^")) (4.75 «) 

depend now on a, and under the normalizing conditions (4.2): $+(00, a) = I and (4.76): Ki(a) < 
1^2 («) < . . . i^k (ct) are determined uniquely. 
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Assume now that for all a from some neighborhood of some point a,, € C the condi- 
tion n G+) = #(7V(^>(a)) n G+) is satisfied (thus the block-entry S""($(a))ii is 
a square matrix for all a which are close to a^. Assume also that det (5*" (<!>(«) )ii) ^ 0, but 
det (S"'($(q;o))ii) = 0. According to Theorem 4.2, applied to the matrix-function $(«) with 
a. 7^ tto) (and a which is close to a„), the matrix function <I>(z, a) {a ^ a„) admits the factoriza- 
tion of the form $(z, a) = ^+{z, a) ^-{z, a), in other words, D{z, a) = I for a ^ aQ. Moreover, 
the factors ^+iz, a), ^-{z, a) depend on a holomorphically for a ^ ao (this may be established 
from the explicit formulas (4.53), (4.54) for the factors). However, according to Theorem 4.1, 
the factorization of this form is impossible for a = ao (otherwise the Lyapunov-Sylvester matrix 
S^{^{a„)) would be invertible). Thus, Z)(ao)) 7^ /, i.e. not all Kj{a(t),j = 1,2, ... ,k, vanish 
(although still X^i<j<fcKj(ao) = 0) . So, the factorization (4.74) undergoes a bifurcation at the 
value of the parameter a such that det S^{^{a)n) vanishes. The natural question arise how to 
describe the bifurcation in a clear way? We believe that for rational matrix function $ in general 
position it may be done more or less explicitly. 

5 . SOME HISTORICAL REMARKS 

There are many different ways to specify and represent analytic functions: for instance, 
Taylor series, decomposition in continuous fractions, representations by Cauchy integrals or 
by Fourier integrals, etc. (Of course, the distinction between different representation methods 
is often artificial and hard to make.) In the first half of 1970's this toolkit was enriched by 
an additional representation method: so called system realizations of analytic functions. The 
sources of the theory of system realizations belong to several different domains, in particular 
synthesis theory of linear electrical networks, the theory of linear control systems, and the theory 
of operator colligations (or nodes) and their characteristic functions. Investigations in these 
theories were carried out by representatives of different scientific disciplines. The investigations 
done by the mathematicians have their root in the pioneering work of M. S. Livsic, who is the 
forefather of the theory of system realizations. 

In the middle of 1940's M. S. Livsic has introduced the notion of the characteristic function 
of a linear operator. This notion was first introduced for nonself adjoint extensions of isometric 
operators with defect indices (1, 1) [LI] and more generally (n,n) [L2, L3] (for n > 1 the charac- 
teristic function is matrix valued), and later for general operators with finite nonhermitian (or 
nonunitary) rank [L4]. M. S. Livsic has discovered the following properties of the characteristic 
function 

1. The characteristic function determines the corresponding operator essentially uniquely up 
to unitary equivalence (first results of this kind are contained already in [LI]). 

2. For each invariant subspace of the operator there is a decomposition of the characteristic 
function into a product of two factors: one of these factors is the characteristic function 

of the restriction of the original operator onto the given invariant subspace, and the other 
factor is the characteristic function of the compression of the original operator onto the 
corresponding coinvariant subspace. 

Property 2 is the so called "multiplication theorem" for characteristic functions. Initially the 
multiplication theorem was established under various additional technical restrictions. The fi- 
nal formulation of the theorem has been obtained in the framework of the theory of operator 
colligations and their characteristic functions. 

^^we formulate these properties, on purpose, in a somewhat rough and therefore sUghtly imprecise form 
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The theory of operator colHgations that was created by M. S. Livsic (and that was further 
developped by his collaborator M. S. Brodskhi, sec the book [Br]) was a natural development of 
the theory of nonself adjoint operators and their characteristic functions. This development was 
also intimately related with applications of the theory of commuting nonself adjoint operators 
to physical problems, in particular to the problems of scattering and to elementary particles 
physics, and later to synthesis problems for electrical networks [L5, L6, LF]. The theory of open 
systems took a definitive shape in the works of M. S. Livsic in the early 1960's; the contents of 
these works were incorporated in the monograph [L9]. 

A theory parallel to the theory of open systems of M. S. Livsic has been developped by 
several other authors under the name of the theory of linear stationary dinamical systems. (For 
an exposition of system theory see [Z], [Fu], [KFA]; we especially recommend the survey [Kaas].) 
The transfer matrix function of such a system is precisely the characteristic function of the 
corresponding operator colligation. As M. S. Livsic has shown in several important examples, 
for scattering systems the characteristic function coincides with the scattering matrix [L5, L6]. 

A different line of investigation leading to the theory of system realizations is connected with 
control theory and with the name of R. Kalman (see [KFA] ) . Here one also encounters the notion 
of the transfer matrix function. Let us emphasize in this connection one important circumstance. 
In physical problems that M. S. Livsic was motivated by there appeared always an "energy 
balance" condition implying J-contractiveness of the characteristic (transfer) function and its 
"symmetry" with respect to the unit circle (or the real axis). It also imposed considering the 
adjoint operator A* together with the state operator A of the system. Energy balance condition 
does not play an important role in control theory and this leads naturally to considering a 
general pair of operators A, B on the state space instead of the pair A, A* . (In fact, R. Kalman 
develops system theory over arbitrary fields rather than over the field of complex numbers.) 

Many results of the analytic theory of electrical networks can be considered as realization 
results for matrix functions of various classes. (There are many expositions of the theory of 
electrical networks; we recommend especially the monograph of V. Belevitch [Bel] that seems 
as if written for a mathematician. See also the survey [EfPo].) The well known Darlington's 
synthesis method for passive networks has been formulated for mathematicians by V. P. Potapov 
in 1966 as a realization problem for passive rational matrix functions [Pot3]; this problem was 
considered in details by Potapov's Ph. D. student E. Melamud [Me]. 

Already in [L4] M. S. Livsic used the theory of characteristic functions for the reduction 
of a nonself adjoint operator to a triangular form, generalizing the theorem of I. Schur that an 
arbitrary matrix can be brought to an upper (or lower) triangular form by unitary equivalence. 
In this approach a multiplicative decomposition of the characteristic function corresponds to 
an "additive" decomposition of the operator itself over a linearly ordered chain of its invariant 
subspaces. This correspondence was used by M. S. Livsic in both directions. Using invariant 
subspaccs of a finite-dimensional approximation of a given operator he decomposes the charac- 
teristic function of the approximating operator into factors, then passing to a limit he obtains 
a multiplicative decomposition of the characteristic function of the given operator, and finally 
using this multiplicative decomposition he constructs a triangular model of the given operator 
which is unitary equivalent to the operator itself. 

A multiplicative decomposition of a meromorphic J-contractive matrix function on the unit 
disk (or on the upper half plane) has been obtained by V. P. Potapov in [Pot3] in a purely 
function theoretic way (the simpler case J = I was handled much earlier in [Potl]). However an 
important special case of V. P. Potapov's theorem has already been obtained by M. S. Livsic in 
[L4] using operator theoretic methods; V. P. Potapov used the investigations of M. S. Livsic as 
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a guideline in his own function theoretic approach. 

It is natural to ask whether one can reduce an operator to a diagonal form by a linear — no 
longer unitary — transformation. In the middle of 1950's M. S. Livsic has posed this problem to 
his then Ph. D. student L. A. Sakhnovich who has obtained numerous results in this direction. 
In particular the problem of the reduction of a nonself adjoint operator to a diagonal form has 
lead L. A. Sakhnovich to a relation of the form 

AX - XA* = GJG* (5.1) 

considered as an equation for X (formula (3) in [SI]). This relation appeared in a hidden form also 
in the theory of J-contractive matrix functions (see formula (22) and the following unnumbered 
formula in Chapter 2, Section 4 of [EfPo]) and in the theory of classical interpolation problems 
(formula (12) in Section 1 of [Kov]). A relation of the form (5.1) in [Kov] appears exactly in 
connection with what we called the chain identity; however the chain identity plays there a 
secondary role and is not highlighted. It should be mentionned that all the main ideas of the 
paper [Kov], published in 1983, have been suggested by V. P. Potapov some 10-12 years earlier; 
unfortunately, V. P. Potapov's contribution is not adequately reflected there. An identity of the 
form (5.1) has also been considered by L. de Branges in the framework of a certain generalized 
moment problem (in the language of Hilbert spaces of entire functions; see Theorem 27 of [Bran] 
and also formula (6.10) in [GolM]). A. A. Nudelman has used an identity of the form (5.1) 
at the basis of an abstract scheme that he developped for considering classical interpolation 
problems. It is interesting to note that an identity of the form (5.1) (and a related identity 
X — AX A* = GJG*) appears in an entirely different context as well, namely in connection with 
fast inversion algorithms for structured matrices (Toeplitz, Hankel, Vandermonde, etc.); see the 
survey [KS]. 

Everything needed for the theory of system realizations in its present form has thus been 
available by the middle of 1970's. An important step was taken by L. A. Sakhnovich in [S2] (a 
detailed exposition of these results is contained in [S3]): he studied the spectral factorization 
of a rational matrix function R with both R and R~^ given as transfer functions of the cor- 
responding linear systems (operator colligations). The spectra of the state space operators A 
and B of these systems do not intersect. One considers a pair of Sylvester-Lyapunov equations 
AT -TB = FiGi, SA- BS = F2G2, where Fi, Gi, F2, G2 are the input and the output 
operators of the systems realizing R, R~^. It is shown that if the corresponding blocks of the 
solutions T and S of these equations arc invertible then the matrix function R admits a spectral 
factorization, and formulas for the spectral factors (having the form of the formulas (4.53-4.54) 
of the present paper) are obtained. N. M. Kostenko (a Ph. D. student of L. A. Sakhnovich) has 
shown in [Kos] that the invertibility of these blocks is necessary for the existence of a spectral 
factorization. Relations which are analogous to zero-pole coTipling relations were considered al- 
ready in [SI] and used in [S2]. Let us notice that the spectral factors constructed in [S2] are 
simply the chacteristic (transfer) matrix functions of subsystems arising by restricting the re- 
spective state space operators to the corresponding spectral subspaces. This is the form taken by 
the multiplication theorem of M. S. Livsic in the current situation (the matrices T and S define 
"metrics" which arc now just bilinear functionals, neither positive definite nor even hcrmitian). 

Unfortunately the paper [S2] did not have the impact it deserved. The subsequent develop- 
ment of the theory of system realizations is connected with the name of I. Gohberg. I. Gohberg 
has also lead and inspired a coherent work of many mathematicians in the theory of system 
realizations and its applications and this theory experienced a fast growth from the late 1970's 
onward. Already in 1979 there appeared the monograph [BGKl] dealing with spectral factor- 
izations of rational matrix functions given as transfer functions of linear systems (operator 
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colligations). The paper [GKLR] considers the realization problem for matrix functions R and 
as transfer functions starting with "local data" (principal parts of Laurent series for R and 
R~^ at each pole); see also [BGRl]. In the beginning it was assumed that the pole sets for R 
and do not intersect; later the general case when these sets may intersect was considered as 
well. These and many other questions are considered in great details in the monograph [BGR2]; 
see also [KRR]. The collection of papers [CoMe] is dedicated to the spectral factorization for 
rational matrix functions based on the theory of system realizations; the papers [BGK2] and 
[BGK3] are especially related with our exposition in Section 4. There are results on system 
realization for rectangular (non-square) rational matrix functions given by local data [BGRa]. 
There are also realization results for matrix functions on a Riemann surface [BV], using deep 
new ideas of M. S. Livsic and his collaborators on characteristic functions for commuting tuples 
of nonself adjoint operators [LKMV]. 

Factorization of matrix functions is a tool for many other problems, e.g., the theory of inverse 
problems for differential equations and prediction theory of stationary stochastic processes. If 
the corresponding matrix function is rational, this factorization (which is a technical tool for 
the original problem) may be carried out using system realizations which then become involved 
in the solution of the original problem as well. See, e.g., [AGl, AG2]. It is clear that the theory 
of system realizations can be successfully used also for the solution of the problems considered 
in [Yag]. It would be interesting to connect the questions considered in [Dei] with the theory of 
system realizations. 
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